





NASA Contractor Report 166061 


OPTIM 

Computer Program to Generate 
a Vertical Profile Which 
Minimizes Aircraft Fuel Burn 
or Direct Operating Cost 

User’s Guide 


Analytical Mechanics Associates, Inc. 
Mountain View, California 94043 


Prepared for 

Langley Research Center 

under Contract NAS 1-1 5497 


NASA 

National Aeronautics and 
Spurn Administration 

Langley Research Center 

Humpton. Virginia .MtWh 


May 


Jiummi “PW-W# 1 ! ■ I ‘ !| ii iA^iana^liNWHIMU r -*msr^ 


FOREWORD 


The development of thin computer program — referred to as 
OPTIM -- was supported under NASA Contract No. NAS1-15497 ,by 
Langley Research Center, Hampton, Virginia. The project technical 
monitors were Samuel A. Morello, Kathy H. Samms, and Robert E. Shanks. 
At AMA, Inc., the project manager was John A. Sorensen, with 
engineering support provided by Mark H. Waters. The project pro- 
grammers were Marianne N. Kidder, Quyen T.L. Nguyen, and Leda C. 
Patmore. 

OPTIM ts an extensive modification of an original program 
developed by Hot nr. Krzbergor and Homer Q. Lee of NASA Ames Research. 
Center. Technical discussions with Dr, Erzberger and Mr. Leo are 
grateful Iv acknowledged. Also, suggestions for program improvement 
bv Ms, Samms and other members of the NASA Langley Research Center 
staff have been greatly appreciated. 

This User's Outdo describes the program input, program out- 
put, and general organization. Appendix A presents the technical 
material upon which the program Is based. Appendix B presents a 
brief explanation of each ot the program subroutines. 
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INTRODUCTION 


This document is a technical description and a user's guide for a 
computer program — called OPTIM — which is used to design near optimum 
vertical profiles for turbojet powered aircraft. Specifically, the program 
generates a profile of altitude, airspeed, and flight path angle as a function 
of range between a given set of origin and destination points for particular 
models of transport aircraft provided by NASA. Inputs to the program include 
the vertical wind profile, the aircraft takeoff weight, the costs of time 
and fuel, certain constraint parameters and control flags. The profile can 
be near optimum in the sense of minimizing: (a) fuel, (b) time, or (c) a 
combination of fuel and time (direct operating cost (DOC)). The user can 
also, as an option, specify the length of time the flight is to span. The 
theory behind the technical details of this program appears in Appendix A. 

OPTIM is an adaptation and an extensive modification of another program 
developed by Dr. Heinz Erzberger and Mr. Homer Lee of NASA Ames Research 
Center for the IBM 360 computer. The present program has been modified 
to operate on the NASA Langley Research Center CDC 6600 and PDP 11/70 
computers, 

A companion program, has been constructed which takes the output of 
OPTIM as input. This companion program — called TRAGEN — simulates an 
aircraft steered to follow the profile generated by OPTIM. The user’s guide 
for TRAGEN appears as a separate document [1] . 

OPTIM has the following applications: 

% 

1. It can determine how much fuel consumption and operating costs can 
be reduced by flying an optimum path rather than a reference tra- 
jectory specified in the pilot’s handbook. 

2. It serves as a benchmark for evaluating sub-optimum algorithms. 

3. It can be incorporated into an airline’s flight planning system. 
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4. It can be incorporated into advanced automatic air traffic control 
software. 

5. It serves as the basis fox. --the -design of an advanced flight manage- 
ment system. 

6. It can be used to assess the advantage of alternate engines or 
aerodynamic changes on air transport operating costs. 

The chapters in this document are organized as follows: 

1. Chapter II explains the meaning of input variables required to run 
the program. 

2. Chapter III explains the meaning of program output variables and 
options. 

3. Chapter IV presents a subroutine layout and flowcharts explaining 
the basic organization of the program. 

4. Appendix A briefly summarizes the theory behind the OPTIM program. 

5. Appendix B explains the purpose of each of OPTIM 1 s subroutines and 
functions. 
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INPUT DESCRIPTION 

To run OPTIM requires up to six ii mt cards and up to three designated 
data files* The meanings of the variables on the input cards are given 
first. The program uses Unit 5 as the card input source. 

Card 1 

This card is the header that appears at the beginning of the run. The 
input has an 8A10 format. 

Card 2 

This card consists of ten integer variables used as flags to control 
the operation of the program. The input numbers are right-justified and 
have a 2014 format. They are: 

ICTAB ICOUT IPRINT I VP I I WIND ICALT ISPUiT I GRAF IAC 

IDCC. 

The meaning of each of these variables is as follows: 


ICTAB In generating the climb and descent portions of the profile, 
the program uses a table of optimum cruise conditions. This 
table, called the "Cruise Table", gives the optimum cruise 
altitude and airspeed as a function of cruise weight. The 
program uses the results of this table to produce boundary 
conditions for the optimization process. The program has 
the option of generating a new cruise table or using a cruise 
table generated by a previous run. A new table should be 
generated either if a different aircraft model (aerodynamic and 
engine data), wind profile, aircraft heading, temperature 
variation, or cost penalty on time or fuel is used. Values of 
ICTAB are: 
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ICTAB = 0: Causes the calculation of a new cruise table. 

ORIGINAL PAGE IS 

OF POOR QUALITY ICTAB - I; Uses an existing cruise table to be read in 

from Unit 8. 

ICTAB = 2: Time-of- arrival fixed (create new cruise tables). 

Creating a cruise table when cruise altitude is free to 
vary is time consuming. (See input option ICALT) . Free 
altitude optimization runs may take up to ten times as long 
as fixed altitude runs because of this calculation. It is 
therefore advisable to save (option ICOUT) and reuse old 
cruise tables whenever possible for free altitude runs 
(ICALT = 0 or 3), It is not possible, however, to save 
or reuse fixed altitude or time-of -arrival cruise tables. 

ICOUT This option controls the writing of a cruise table on Uni;: 8, 

It is only operative during free cruise altitude runs. The 
user is responsible for saving Unit 8 output in a permanent 
file after completion of a run, for rewinding it between 
runs, and for recovering before a new run. Values of ICOUT 
are; 

ICOUT =0; Do not write cruise table on Unit 8. 

ICOUT =1; Do write cruise table on Unit 8. 

1PRINT This flag is used to control the amount of printout during 
the program computation process. Values of IPRINT are: 

IPRINT = 0: Normal mode (see Chapter III, for output description) 

IPRINT = 1: Extra printout included. This produces detailed 
output useful for debugging cruise table 
calculations. Some familiarity with thr program 
is necessary in order to use this output. 

IPRINT ® 2: Less printout (input, cruise table summary, summary \ 

of flight parameters (time, distance, fuel) for 
each climb and descent iteration, optimum profile 
summary) . 

IPRINT = 3: Minimum printout (input, cruise table summary, 

summary of flight parameters for optimum profile.) 


IVPI 


IWIND 


1C ALT 
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To determine the optimum profiles, the user has the option 
of using airspeed as the control variable (with thrust fixed) 
or both airspeed and thrust as control variables. If thrust 
is fixed, it is set to the maximum value for climb and to the 
minimum value (idle throttle) during descent. Using two con- 
trols gives slightly better performance trajectories in terms 
of lower overall cost. However, for short-haul operations, 
the profile shape can be substantially different, as shown in 
Appendix A. Values of IVPI are: 

IVPI = 0: Optimize using only airspeed as a control. 

IVPI « 1: Optimize using both airspeed and thrust as 

controls. 


An arbitrary wind profile can be read in on Unit 7. It 
gives the wind speed and heading as a function of altitude 
(above sea level). Values of IWIND are: 

IWIND =0: No wind used. 

IWIND = 1: Input wind profile varies with altitude but 

is constant over the entire range of flight. 
(See description of input wind data, p. 9 ). 

IWIND = 2: Separate wind profiles are input for climb, 

cruise and descent portions of the flight. 


The program has the option to generate a three-part profile 
(consisting of climb-cruise-descent) or a two-part profile 
(consisting of cruise- descent) with or without fixed cruise 
altitude. An optioru-tQ add .a_step climb segment during 
cruise is also available. The values of ICALT arc: 


ICALT = 0: 
ICALT * 1 : 


ICALT = 2: 


Three-part profile with a frer cruise altitude. 

Three-part profile with fixed cruise altitude. 
Wi,th this option, if the input range is not long 
enough to allow the aircraft to climb to and 
descend from the input cruise altitude, a 
feasible altitude will be sought, ne final 
altitude will be some multiple of 2000 feet less 
than the input altitude. 

Three-part profile with fixed cruise altitude 
and step climb. The program will assume a 
4000 foot climb at maximum thrust after attain- 
ment of the fixed-cruise altitude. The optimum 
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distance into cruise at which the step climb 
starts is solved for along with the other 
optimization variables. 


ICALT * 3: Two-part profile with free cruise altitude. 

Initial cruise weight and range-to-go are 
input. OPTIM solves tor the optimum initial 
cruise altitude and airspeed in addition to 
the rest of the profile. 


ICALT = 4s Two-part profile with fixed cruise altitude. 

Initial cruise weight, range-to-go, and 
altitude are input. OPTIM solves for the 
optimum cruise airspeed with altitude fixed. 

ICALT ■ 5: Two-part profile with fixed cruise altitude 

and step climb. This option is similar to 
ICALT 9 2, except that the flight starts at 
an input initial cruise weight and a ltitude. 


ISPIMT This flag allows the user to remove the 250 kt indicated 
airspeed limit below 10,000 ft. Values of ISPLMT are: 


ISPLMT =0: No V IAS limit for altitude below 10,000 ft. 

ISPLMT - 1; 250 kt V^g limit below 10,000 ft (nominal). 

ISPLMT = 2: 250 kt V limit below 10,000 ft for descent only. 

IAd 

IGRAF This flag controls the output of a data set containing the 
optimum profile which is used for generating graphs and 
as input to the TRAGEN program. The data set is output 
on UNIT lie Values of IGRAF are: 


IGRAF =0: Do not write an output data set. 


IGRAF = 1: Write a data set using Unit 11. 

IGRAF * 2: Write a data set using Unit 11. Printer plots 

are drawn for the four variables: Mach no., 

flight path angle, altitude and fuel versus range. 


IGRAF >2: A data set is written on unit 11. Printer 

plots are generated for all variables. 


IAC This flag is used to select which aircraft model to use to 

generate the optimum flight profile. Current values of IAC 
are; 
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I AC - 2; Medium-range three-engine jet transport air- 
craft. 

IAC 3: Medium-rant a two-engine jet transpo rt aircraft. 

Descent may be constrained to a constant sink rate to main- 
tain cabin pressurization at a given differential until cabin 
pressure reaches sea level pressure. At present this differ- 
ential is set at 10 psi (pounds per square inch). The con- 
stant sink rate is 500 ft/min. Values of IDCC are: 

1DCC - Os Unconstrained descent 

IDCC = 1: Constrained descent. 

Card 3 (Optional) 

This card has four real variables with format 8F10.2. It is used only 
when a new cruise table is to be generated (ICTAB * 0 or 2 on Card 1) • ihe 
variables are: 

FC TC DTEMPK PSIG. 


c r ; ^ 4 - 
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IDCC 
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The meanings of these variables are as follows: 


FC 

TC 


DTEMPK 


PSIG 


This is the cost of jet fuel in $/lb (e.g. , 0.15). 


This is the cost of time in $/hr (e ^_. t _,j>QP.00) . Both 
FC and TC are used in the cost fur tion which is minimized 
by the program. If t ime-of-arrival is fixed, this variable 
is ignored. 


This is the temperature variation from standard atmospheric 
conditions in degrees Kelvin. 

This is the aircraft ground heading in degrees. It is used 
along with the wind heading to compute aircraft heading with 
respect to the airmass. 


Card 4 (Optional) 

This card has three real variables with format 8F10.2. It Is used 
only when a new cruise table is to be generated ( 1CTAB = 0 or 2 on Card 1) 
The variables are: 
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WEIGHT WN DEW. 
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The meanings of these variables are as follows; 

WEIGHT This is the maximum value that the weight of the aircraft 

can be in pounds (e.g. , 150,000 lb). The first cruise table 
will be generated at this weight. 

WN This is the minimum value that the weight of the aircraft 

can be in pounds (e.g., 110,000 lb). 

DEW This is the incremental cruise weight in pounds between each 

table (e.g., 5,000 lb). Starting with WEIGHT, a cruise table 
will be generated for cruise weights of WEIGHT, WEIGHT-DEW, 
WEIGHT-2 DEW, etc., down to WN. A maximum of ten cruise 
tables can be generated because of internal program array 
dimensioning. 

Card 5 

This card has five real variables with format 8F10.2. The variables 

are; 

WTO RANGE DEIN HCRUZ TEND. 

The meanings of these variables are as follows: 

WTO Aircraft takeoff weight for three-part profile or initial 

weight for two-part profile in pounds (e.g., 136,000 lb). 

NOTE: This value must be less than or equal to WEIGHT of Card 4. 

RANGE Rarige-to-go in nautical miles (e.g., 200 n.mi.). 

DEIN Incremental specific energy in feet between points on the 

optimum profile. As explained later in Appendix A, the 
program uses an energy state method to generate the optimum 
trajectory. Specific energy is the independent variable, 
and varying the size of DEIN affects the smoothness and 
accuracy of the generated profile. If DEIN is input as 0. , 
it is set equal to a nominal value of 500 feet. 
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This is the value of the fixed-cruise altitude (in feet) 
that is used when ICALT is set other than 0 or 3. 

This is the value of the desired time- of -arrival (in seconds) 
that is used when ICTAB is ^et to 2. 


This card has four real variables with format 8F10.2. The variables 


Initial aircraft altitude in feet. 


VTO Initial indicated airspeed in knots. 

HOLNDG Final aircraft altitude in feet. 
VOLNDG Final indicated airspeed in knots. 


Oft ^ 


Note that if a cruise table already exists, only Cards 1, 2, 5, and 6 are 
required. An example of cards 2 through 5 is: 


0 0 3 0 0 

0*15 <600* 

100000* 70000. 

100000 . 1000 . 

0 . 210 . 


1 0 3 

0. 90. 

5000 „ 

500 ♦ ^33000. 

0 . 210 . 


5000. 


In addition to the card input, there are up to three data sets that 
may be used by the program. These are: 


Unit 7 - Wind Data (Optional) 

This data set is used when IWIND is set to 1 or 2. The input consists 
of the magnitude of the wind and the direction of its source as a 
function of altitude. The data format — is — . . . , ^ 

If IWIND - 1, a single wind profile applicable to the entire flight is 
read in. This profile consists of a set of n cards. Each card has four 
variables . 

HWIND(I) PSIW(I) VW(I) IE 

There is one card for each I«l,2 N, where N is the number of al titudes used 

for a given wind profile. The meanings of 


9 



original pass m 

OF POOR QUALITY 


Beginning (lowest) altitude at which direction PSIW(I) 
and magnitude VW(I) apply. The program will interpolate 
for values of PSIW and VW when using altitudes between 
HWIND(I) and HWIND(I+1) . 

Direction of the wind vector source in degrees (i.e., 270° 
represents a wind from the West). 

Magnitude of the wind vector in knots. 

End-of-wind- table indicator. If IE * 0 the program will 
expect to read further wind data. If IE - 1 the program 
assumes that a complete wind table has been read in. Note 
that when IE * 1, the corresponding altitude should be equal 
to or greater than any altitude the aircraft is expected to 
reach. 

If IWIND = 2, three wind profiles are read in, one each for climb, cruise, 
and descent (in that order). Each profile is as described under IWIND * 1. 
Each profile must end with a non-zero value for IE. 

Unit 8 - Cruise Table Data 

Data are generated by the program and written on Unit 8 when ICOUT is 
set to 1. Most of this data are also output to Unit 6 when IPRINT « 0. 

Data generated during a previous run are read in from Unit 8 when ICTAB is 
set to I; in this case, only the summary of the cruise table is output to 
Unit 6. The cruise table data that are written on Unit 8 are as follows: 
(the Format is 8E15.7 unless specified otherwise): 


HWIND(I) 

PSIW(I) 

VW(I) 

IE 


Line 1: FC, TC, DTEMPK, PSIA obtained from Card 2. 


Line 2: 


Line 3: 


WUSE - aircraft cruise weight obtained in going from WEIGHT 
to WN in steps of DEW as defined by Card 3. 

HALT altitude - 10,000 to 40,000 ft (or to H mav ) in 
steps of 1,000 ft. 

QiAKLAS minimum drag airspeed - kt 
FBIAS maximum cruise airspeed - kt 

OPMIAS indicated airspeed - kt 


Optimum cruise con- 
ditions at altitude 
H for minimum DOC 
per n.mi. 


ORICWal v 

op pcoh c ‘;;, 


OPMTAS true airspeed « kt 
OPMACH Mach 
EPRS EPR setting 
FDTOPT cost in $/n.mi, (A) 

FUELDt fuel flow rate for FDTOPT 


Optimum cruise con- 
ditions at altitude 
H for minimum DOC 
per n.mi. 


These lines are repeated for altitude varying from H 

min 


Line 4 : ENDATA 


Line 5: 


OPTMAK 
OPT IAS 
OPTTAS 
EMCOST 


FUELD1 


-1G , code for end of a weight set. 

optimum altitude where A (FDTOPT from line 3) 
is minimum 

optimum Mach (at HOPT) 
optimum Indicated airspeed - kt 
optimum true airspeed - kt 
minimum - $/n.rai. 
optimum EPR setting (at HOPT) 
optimum cruise energy (at HOPT) 
optimum fuel flow rate. 


An example of lines 1-5 appears as an output table in Chapter HI. 
Lines 2-5 are repeated for each weight, from WEIGHT to WN in steps of DEW. 


Line 6: 


Line 8; 


Line 9; 


ENDATA 


Line 7 : IWMAX 


-10 , code for end of weight tables. 

(Format 14) number of cruise tables 
(~ ( WEIGHT- WN) /DEW) H.) 


(DLLDEL ( 1 , J ), DLLDLE( 2, J) , J— 1, IWMAX) — the coefficients 
of the derivative of the optimum cost EMCOST as a function 
of energy for each weight table. 

WS (I) cruise weight - lb \ 

EOPTS(l) optimum cruise energy for WS(1) - ft I 

EMSTAR(I) optimum Mach I 

HSTARS(l) optimum altitude - ft | 

PISTRS(I) optimum El'R V Repeat 

/ 1 - 1 . 

ELAMBS(l) minimum cost - $/n.mi. 1 

VTASOP(I) optimum true airspeed - kt 1 

FUELFUD fuel flow rate at PlSTRS(l) - lb/hr. I 


Repeated for 
I - 1, IWMAX 
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OUTPUT DESCRIPTION OF POOR QUALITY 

The output of OPTIM is lengthy, in tabular form, and generally self- 
explanatory. This section presents one sample of each of several different- 
types of output written to Unit 6 • The »ju«.n u ity of this output is controlled 
by the input flag IPRINT. Unit 8 (input flag ICOUT) and Unit 11 (flag 
IGRAF) can be used to store data sets for later use if desired. 

A summary of the flags and input variables is printed at the beginning 
of each run. Table 1 is an example of this output. 

The next table output is the vertical wind profile. This is printed 
when there is a non-zero wind, and the flag IWIND is set to 1 or 2. An 
example is presented here as Table 2. 

The next output tables are the cruise tables. These are printed if 
they are generated as part of the run. That is, if the flag ICTAB is set 
to 0 or 2, new cruise tables are computed based on other input data. 

Table 3a is an example of a cruise table for a cruise weight of 100000 lb. 
Each column of this table is as defined on pp. 10-11 (in the explanation 
of Line 3, Unit 8). The format of the printed version varies slightly 
from the format of the data set output to Unit 8. 

If the cruise tables are not generated as part of the run, but read 
in from Unit 8, then a table such as Table 3b is printed. This table pre- 
sents the derivative of the cost per n.mi. (referred to as X ) a a func- 
tion of change in the cruise energy. The coefficients of this derivative 
are shown for cruise weights from WN to WEIGHT in steps of DEW. (See input 
description, page 8). 

Following Tables -3a or 3b is a summary of the optimum cruise condi- 
tions for each of the distinct cruise weights in terms of minimum direct 
operating costs. This table is obtained by interpolating the conditions 
of each cruise table, such as Table 3a, tc determine the optimum altitude/ 
airspeed/power setting combination for a given cruise weight. An example 
of this summary table is shown as Table 4. 
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Table 1. Summary of Input Flags and Data 
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Table 2. 


^ 4 ~ i.7 

OF POOR QUALITY 

Vertical Wind Magnitude and Direction as a Function 
of Altitude. P31W is direction of wind source. 


OWIND DATA CLIMB 
ALT (FT > » VW ( KNOTS ) » 
0 . 0.00 
40000. 25.00 

OWIND DATA CRUISE 
ALT (FT > » VW ( KNOTS ) » 
0 . 0.00 
40000. 50.00 

OWIND DATA DESCEND 
ALT (FT ) » VW ( KNOTS > » 
0 . 0.00 
40000. 100.00 

0 AIRCRAFT HEADING » 


VWCFT/SEC) » 
0.00 
42.20 

VW (FT/SEC ) » 
0.00 
84.3? 


PS1 W < DE G> 
90. 

? 0 , 

PS I W ( DEC > 
90. 

90. 


VW (FT/SEC ) t 
0.00 
166.78 
90. DEG 


PSIW(DEG) 

90. 

90. 


After the summary table is printed, an estimate is made of the initial 
cruise weight (takeoff weight minus fuel burned during cximb). This value 
is then used to generate a new line of the cruise table based on inter- 
polation. The result is printed out and is shown here as Table 5. The 
variable LAMBDA from Table 5 is used as a key search variable in constructing 
the optimum trajectory. This is explained in more detail in Appendix A. 

After the first five tables are printed, the program goes into an 
interative search process to compute the optimum climb and descent portions 
of the trajectory. The number of iterations varies and depends on the range 
to be flown, whether a two or three-part trajectory is solved for, whether 
optimization uses airspeed or airspeed/thrust as controls, and how close 
to the desired range the final trajectory is supposed to be for conver- 
gence, An example of the output given in tabular form for the climb 
trajectory is shown in Table 6a. 

The descent trajectory is computed backwards in time starting with 
the estimated landing weight and ending with the final cruise conditions. 

An example of this descent trajectory is shown as Table 6b. 

A summary of the climb, cruise, and descent segments of a three— part 
trajectory in terms of fuel, distance, time, cost and cost/n.mi. is printed 
finally in the form of Table 7a. This table is the essence of the 
program's output. It shows how close to the desired range the trajectory 
came, and what total cost resulted. This tabic is used to rapidly compare 
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Table 3b. Derivative of the Cruise Cost with Respect to 
Cruise Energy. (IPR^NT - 0 or 1 only) 
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Table 4. Cruise Table Summary. 
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Table S . Interpolation of Cruise lable for Estimated Weight 
at Top of Climb. I PRINT - K 
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Table 7a. Summary of Optimum Results for a Climb-Cruise- 
Descent Profile. 
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131262. 

tas 

482. 

482. 

COST! t /NHI 
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Table 7b. Summary of Optimum Results for a Cruise-Descent 
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one run against another. An example of a two-part trajectory (cruise- 
descent) is shown in Table 7b. Tables 7 are output for any value of 
the printout flag IPRINT. 

Table 8 is an example of the table of optimum trajectory variables 
which may also be written to a data set for use in additional plotting 
routines or for the program TRAGEN if IGRAF >_ 1. Note that this summary 
shows GAMMA, which is the flight path angle with respect to the air mass, 
rather than the FPTH shown in Tables 6a and 6b. FPTH is the flight 
path angle with respect to the ground. 

Table 9 is a summary of the cruise performance which is printed at 
the end of each run in which a cruise section is produced. Table 9 
shows the steady cruise-climb conditions every 100 n.tni. to cover the 
expected range of the flight. Note that as fuel is burned off, the 
optimum altitude rises, and the optimum power setting and cost per n.mi. 
change. If cruise altitude is fixed, this will be indicated in the 
various cruise tables. Table. 9 may also be written to Unit 11. 

Table 10 is an example of the step climb portion of the profile 
calculated when step climb is included in the optimization. 

Figures la- lb are examples of plots produced when IGRAF > 1. 


Unit 11 Output data set (optional) 

The variables which describe the optimum vertical profile followed 
by the aircraft are written as output on Unit 11 when IGRAF » 1. Printer 
plots are also generated when IGRAF is greater than 1. 


The output data are used for two purposes: (a) they serve as in- 
put to plotting routines so that a more convenient record of the data 
than tabular listings can be obtained, and (b) they serve as input of 
points on the nominal optimum trajectory for the program TRAGEN described 
in Ref. 1. 
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Table 10 


Step Climb Portion of the Optimum Profile. (IPRINT < 3) 
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,729 

414. 7W 

1.941 

724.41? 

7.237 

414.412 

114.415 

11.971 

14.414 

*3943.49 

1144.911 

41444.449 

14444.444 

. 774 

419.741 

1.49? 

774. 129 

7.23* 

419.547 

144.447 

14.944 

19,174 

*37*3.14 

1979.179 

44129.124 

14544.444 

.724 

413.273 

t.7ll 

323. *14 

7.737 

419.11* 

174.554 

74.44? 

13.934 

*3737.9 7 

*49.19* 

44423.149 

17444.444 

.774 

415.137 

1 .141 

344,173 

7.737 

419.424 

212.334 

74.410 

12.434 

*3474.41 

939.4*1 


The output has up to six binary records of the form: 

Record 1; WORD, NWORD 

WORD may be: CLIMB, CRUISE, OR DESCEND, 

NWORD is the number of points stored for the 
specified flight segment. 

Record 2: An NWORD by 12 matrix. For example, for Climb, 

Record 2 contains for JCL1MB * 1,.., NWORD: 


CGRAF ( J CLIMB ,1) 

= 

E 

CGRAF( JCLIMB , 2) 

* 

J 

CGRAF(JCLIMB,3) 

* 

MACH 

CGRAF( JCLIMB, 4) 

s 

VTASK 

CGRAF( JCLIMB, 5) 

= 

GAMMA 

CGRAF (JCLIMB, 6) 

s 

FUELUZ 

CGRAF (JCLIMB, 7) 

- 

EPR 

CGRAF (JCLIMB, 8) 

= 

0 

CGRAF( JCLIMB, 9) 

= 

TIME 

CGRAF (JCLIMB, 10) 

= 

DIST 

CGRAF( JCLIMB, 11) 

= 

HDOT 

CGRAF (JCLIMB, 12) 


VIASK 


Specific energy - ft 
altitude - ft 

Mach 

true airspeed - kt 
flight path angle - deg 
fuel burned - Eb 
EPR setting 
blank 

time - sec 

range traveled « n. mi. 
vertical rate - ft/min 


indicated airspeed - kt 


Flu? same variables are stored in WRAP (JDESCN,J) for the descent portion, 
and SCRAP for the combined cruise and step climb portions. 
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Figure la. Printer Plot Example. Altitude vs Range 



Figure lb. Printer Plot Example. Mach Number vs Range 


IV 


PROGRAM ORGANIZATION 

This section gives a brief overview of the process used in QPTIM to 
generate an optimum vertical profile for aircraft operating between two 
fixed points. The overview is in the form of steps and program flow 
charts. Then a brief description is given of each of OPTIM 1 s 62 subrou- 
tines and 19 functions. Charts showing the interrelations between these 
subroutines are also presented. 

The technical details upon which the program is based are presented 
in Appendix A. A more detailed description of the subroutines is presented 
in Appendix B. More details concerning the program’s origin and concepts 
upon which it is based can be found in Refs. 2-6. 

OPTIM has been configured with a short main program which calls the 
input subroutine, ALLIN, and the major control subroutine, OPTM56. 

OPTM56 then follows one of two paths depending on the input parameters. 

It may synthesize a fixed-range, two-or three-part trajectory. Or, it 
may call OPTTOA to determine a trajectory meeting a fixed time-of-arrival 
constraint. 

The fixed time-of-arrival is accomplished by iterating on the cost 
of time TC between passes of the program. In other words, an outer loop 
is used to iterate on the cost-of-time coefficient TC so that a fixed 
time-of-arrival is achieved. The inner loop determines the optimum ver- 
tical profile with a f ixed-cost-of-time set by the outer loop and a 
fixed cost-of-fuel FC which is input. 


Figure 2 shows a flow chart of the steps followed by OPTIM to synthe- 
size a fixed-range, two- or three-part trajectory consisting of climb, 
cruise, and descent profiles. For a case where f ixed-time-of-arrival is 
desired (ICTAB=2), Fig. 2 represents the inner loop of the program. (The 
outer loop is discussed later.) For the input flag ICTAB set at 0 or 1, 
this is the normal program flow. The steps followed are: 
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Read all input and place it in COMMONs/INPUT/ and /CRUISE/. 

(The latter is filled by input only if an old cruise table is 
being read.) 

Generate the cruise table if it has not been read in. 

Test on trajectory type (input variable ICALT). If this is 
a two-part trajectory, go to (5a). 

Set P . s 1.00 or 0.0 based on flags IVPI and ICTAB. Use this 
min 

and other quantities to generate an optimum trajectory of range 
R . If R is less than R , go to (5). If R is within c of R , 
then the desired trajectory has been generated. If R is greater 
than R max > then the trajectory achieves optimum cruise altitude. 

If so, compute the cruise distance d c = R - (d^ + d^) , and 
compute the final cruise weight. Next, use this updated cruise 
weight to recompute a refined descent trajectory. Use this 
refinement to complete the three-part trajectory. 

For a two-part trajectory, use A ^ and input quantities to 
generate the optimum descent profile. Compute the cruise 
distance d c - R - d^, and compute the final cruise weight. Use 
this updated cruise weight to recompute a refined descent trajec- 
tory. Use this refinement to recompute the two-part trajectory. 
(Note, if initial altitude is not specified, the starting optimum 

cruise altitude and X are chosen based on the input initial 

opt 

cruise weight. If initial cruise altitude is specified, this is 
fixed, and A t is obtained from the cruise table corresponding 
to the initial altitude and weight.) 

Compute P that causes cruise altitude just above 10000 ft. 
r max 

(See Appendix A) , Use the input quantities and P to generate 

max 

the optimum trajectory of ran^e R ^ • Compare inpi t range R with 
R ^ ^ ^ * S ^ ess t ^ ian then no trajectory is computed. 

If R is within e of R m i n > then the desired trajectory has been 
synthesized ( t - 5 n.mi.). If R is greater than R ^ , then go 
on to (6). 


For R 


< R < 


R 

max 


use (P 


R . ) and 
min' 


min max' ' max 

starting points to compute subsequent values of P 


(F . , R ) as the 
min’ max 


Iterate 


on P until the desired range is achieved. (Again, refer to 
Appendix A for more details.) 
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In steps (2), (3), and (4) of Fig. 2, an optimum trajectpry is 
generated based on the parameter P. Computing the points on 
an optimum profile and evaluating the results for a fixed value of X 
consists of nine steps which are presented in flow chart form in Fig, 3. 
The input quantities are P # the initial weight W^, and the initial 
and final values of specific energy (E^ t E^) (which are computed from 
initial and final altitude and airspeed (h^, , (h f , V^)). The steps 

follow the analytical expressions developed in Appendix A, and they are: 


1, The climb iV** 1 hu. ned is estimated based on an empirical 

equation whic. < ■ -fiction of initial and final energy, the cost 

parameters ci and tiiu* , C t ) , and the initial weight. 

The climb fuci p stim^e^ls^aubtracted from the initial weight 

to obtain an estxmw.e oJL-the initial cruise weight W ci . 

Based on W the optimum cruise cost X is obtained from 
ci opt 

the cruise table. This X term is multiplied by a percentage P 

opt 

to obtain the cost parameter X^ used for the optimization of the 
climb profile. This modified cost term and W ci are used to 
interpolate in the cruise table to obtain the estimate 
of the initial cruise energy E^. If a time -of -arrival or cruise 
altitude are fixed, then P = 0, and X - X Qr%t * 


2 . 


The optimum climb trajectory is generated by stepping along 
in energy from E^ to E^, l* 1 steps of AE. At each point, the 
Hamiltonian 


up 


Min 

V.TT 


C r f + CL - X 


i 


)■ 


E > 0 


is minimized by choosing the best value of airspeed v (TVTI = 0) 
or airspeed and EPH ir(TVri = 1). Constraints limit, the 
range of airspeed and thrust over which the search takes 
place. Energy and altitude are monotonically increasing. 

Based on the change in energy and airspeed, the approximate 
changes in altitude, flight path angle time t, fuel 
burned F, and distance traveled X are computed. This pro- 
cess is stopped when energy E ^ is reached, if the weights at top 
of climb is more than 200 lb off the estimate W p this step is 

represented with improved values of X^. 
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Compute P for h - 10000. 

IRAK MX 

Synthesize trajectory for specified 

%.*• V h - v - h ** v * 


Profile synthesized 


'max* "V h f V "f ’f 
Range of synthesized profile -R * 
Store P ■ into the P (*) array" 1 n 
m Into the R (•) array 


m 

^ R| R , 



min 





o profile existed 


Use(P , R in ) and (P R *) as starting points 
v max min min max 

Compute P. Synthesize profile for specified 

P, W., h 1 , V t , h ff V f 

Range of synthesized profile - d 

Store P Into the P (•) array and d into the R (•) array 


Select three best 
points from the P (0 
and R (*) arrays 


Profile synthesized 


g'l— -■ Cruise 
iMtitude 
Fixed 
X ? / 


Yes 


Choose X , h, 
opt i 

as function of W 


frorr, Cruise Table 


Choose X from Cruise Table 
c7ion~o h 


Generate descent profile 
Compute cruise distance and 
final cruise weight 
Refine descent profile 



Figure 2. Continued. 
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Figure continued 
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3. The final cruise weight is next estimated. For the 

climb-descent type of trajectory (IVP1 « 1; see Appendix A) no fuel 
is burned in cruise , is set equal to W^, and the program 
proceeds to Step (4). Otherwise, the program uses empirical 
equations to compute cruise distance d^. Then, with 

initial cruise cost X^, fuel bum rate, and estimated ground 

speed Vg, it computes an initial estimated fuel burned during 

cruise. This value is used to estimate the aircraft weight 

W half-way through cruise. A new cruise cost X and fuel 
c * 

burn rate f are obtained from the cruise table. These values 

are used to refine the estimate F of the fuel burned and the 

c 

final cruise weight W cf . 

4. The estimate W cf is used to obtain from the cruise table the 

final optimum cruise values of cost X £ . The value of X,r 

ct f 

used for descent optimization is found by multiplying X^ by the 
percentage P. The final value of cruise energy used for descent 
optimization is found by interpolating from the cruise tables 
using W £ and A f . 

5* The landing weight is estimated using empirical equations. 


6. Similar to Step (2) , the optimim descent trajectory is 
generated by stepping along in energy (backwards in time) 

from to E^, in steps of AE. At each point, the Hamil- . 

tonlan 



Min 

XJL 


c f f + C t - Xf V, 

fET^ 


E < 0 , 


is minimized. Again, changes in altitude, flight path angle, 
time, fuel burned, and distance traveled during descent are 
computed. 

7. The cruise distance is refined. For aircraft reaching 

optimum cruise altitude (P - 0) , cruise distance d is R - 

c 
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d - d, . Otherwise, a formula using the final values of 
up an 

H , H , and the slope di^/d£ are used to compute d . With 
up an 

this more accurate distance, improved values of and W f 
are computed. 

8. The descent profile is recomputed based on the improved value 

of W^. This produces a better value of descent distance d^ 

which in turn produces a better value of cruise distance d c - 

9. Based on Steps (l)-(8), a table (Table 7a) is generated which 
summarizes the distances traveled, end conditions, fuel burned, 
and costs of each segment of the trajectory plus the trajectory 
as a whole. 

The above steps and flow diagrams of Figs. 2 and 3 are brief summaries 
of the process taken by the OPTIM program to implement the vertical pro- 
file optimization techniques outlined in Appendix A. The reader who is 
interested in more program details is referred to Ref. 3. Further under- 
standing will come from use of the program and study of the individual 

subroutines. 


Two other variations to the basic process of generating optimum vertical 
profiles exist within the structure of OPTIM. The first is the ability to 
simulate a step climb during cruise. This climb is currently 4000 ft from 
one fixed cruise altitude to another (e.g., 33000 ft to 37000 ft). For this 
option, the step climb is computed based on use of maximum climb thrust, and 
the speed is ramped from the optimum at the lower altitude, to the optimum 
at the higher altitude. The climb computation consists of eight 500 ft steps, 
with the aircraft trimmed so that flight path angle rate for each step is zero. 
Then, optimum cruise and descent are computed from the higher cruise altitude. 
The program is set up to compute the optimum-point to begin the step climb if 
this option is chosen. 

The second vari.it Ion is t o constrain the rate of descent at the top of the 
descent portion of the profile. This option might be used to account for cabin 
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pressurization constraints. Currently, OPTIM has the ability to constrain 
the rate at 500 ft/min down to 28000 ft. Then, an optimum descent profile is 
computed from this point downward, where X^ is based on the aircraft weight 
at 28000 ft. 


The capability to generate a vertical optimum profile with fixed range 
and fixed time-of-arrival is governed by the outer loop logic shown in Fig. 4. 
The program begins (as also shown in Fig. 2) by reading in control flags and 
other trajectory characteristics data. If the control flag ICTAB is set to 0 
or 1, OPTIK generates a minimum cost optimum profile based on the input eost-of- 
fuel FC and cost-of-tirae TC. This is the normal mode of operation. 


If the flag ICTAB is set to 2, the program will iterate on the value 
of cost of time TC until the time-of-arrival is within some tolerance_ 
of the desired time-of-arrival T ^ ^ (input TEND) . The logic for this itera- 
tion scheme is shown as Blocks 2, 3, and 4 in Fig. 4. 


The first step is to set TC to zero and generate an optimum vertical 

profile. This profile will correspond to a minimum fuel path. The logic 

to generate this profile is essentially the same as was indicated in Figs. 2 

and 3. The final time T f is recorded. If T fo is greater than T end the 

initial flight profile took too much time. Therefore, time must be penalized 

with positive cost TC. The program then uses logic indicated in Block 3a. 

If T r is less than T the initial flight profile was too fast. Then, 
fo ena 

time must be penalized with negative TC. The program then uses logic indi- 
cated in Block 3b. 

For positive TC, the program generates two more optimum profiles cor- 
responding to costs of time of $300/hr and $600/hr. (If $600/hr produces 
too slow a profile, the program uses $900/hr, etc.) From each of these, the 
times-of-arrival T fl and T f2 are recorded. Then, the program proceeds to 
Block 4 to solve for a TC which will yield the desired T^. 


For negative TC, the program next solves for the minimum cruise value 


of TC, or 


TC 


min 


-f 


rain 


FC. 
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FC is Che input cost of fuel, and f m<n is the minimum cruise fuel rate 
at a particular cruise altitude. The aircraft should not go slower in 
cruise than the cruise speed corresponding to ^ m £ n * Thus, repre- 

sents the boundary on negative values of TC. For this value of TC, the 
variable X Qpt (see Fig. 1 and Appendix A) is zero. A minimum fuel rate 
profile is then generated using T c m ^ n » and the corresponding time of 
flight Tf^ is recorded. (If Tf^ is less than T ^ the difference 
should be made up in a holding pattern at the end of cruise.) For 
greater than T the next value of TC is set by linear interpolation 

between 0 and TC to atteinpt to produce a equal to 

If equals T en< ji the program is finished. If not, Tf£ is 

recorded, and the program proceeds to Block 4. 

In Block 4, three values of time-of-arrival T f (Tf Q > Tf^, Tf 2 > are 
used with three corresponding values of cost-of-time TC (TC^, TC 2> TC^) 
to form the quadratic relationship 

T f « a(TC) 2 + b(TC) + C . 

This equation is solved for TC so that Tf "T^^. Then^the program is 
rerun with this new value of TC used to generate the optimum profile. 

The resulting Tf is compared to T e ^. If it is within e (10 sec), 
the program stops. If not, the new values of (TC, Tf) are used with pre- 
vious values to recompute the quadratic relationship, and the process is 
repeated. This continues until the generated flight paths converge to the 
desired time-of-arrival. 

OPTIM is programmed in FORTRAN, and it consists of the main executive 
program, sixty-two subroutines, nineteen functions, and three block data 
routines. These subroutines and functions are called to execute the steps 
depicted in Figs. 2-4. Brief explanations of the program and its subroutimes= 
are presented in Appendix 13. 

The program subroutines and functions can be grouped into four cate- 
gories: 
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1. trajectory optimization and generation, 

2. aerodynamic and propulsion models, ORIGINAL PAGE fS 

OF POOR QUALITY 

3. flight condition models, and 

4. utility programs. 

Under Category 1, the routines are: 

OPTM56 Serves as the driver program after all input i9 in. 

OPTTOA Serves as the outer-loop driver program when searching on 

TO to achieve the desired time-of-arrival. 

ALLIN Reads ail input. 

CLIMB Generates and stores the climb profile. 

CRUISR Computes a cruise segment from the lower cruise table before 
a step climb. 

CRUISX Compute a cruise segment from a given starting weight proceeding 
for a given range, 

CRZ0P5 Generates the cruise table. 


CTABLE Interpolates in the cruise table for a set of parameters 
at a given weight. 


DESCND Generates and stores the descent profile! 

DESPC Modifies the find portion of the cruise and descent profile 

to provide for constrained descent. 

DRAGC Computes the drag force. 

ESTCD Calls appropriate routine to estimate the cruise range. 

ESTCD2 Estimates cruise range for the tri-jet as a function of 
percent P of \ 

opt 

ESTCD3 Estimates cruise range for the twin-jet as a function of 
cost of time, cost of fuel, and P. 

ESTDF Call appropriate routine to estimate the descent fuel. 

ESTDF2 Estimates the descent fuel for the tri-let as a function of P. 

ESTDF3 Estimates the descent fuel for the twin-jet as a function of 

cruise and landing energy, cost of time, and cost of fuel. 

ESTEP Computes the value of the next energy step during climb 

and descent. 
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FBOUND 


FCLIMB 

FCLMB6 

FCOST 

FDSCN6 

FDRAG 

FOPT 

FTHRST 

FULEST 

FULEST2 

FULST3 

PCCMP5 

PILIMT 

PRETBL 

PRFTBL 

PROFIL 

PRSUM 

STEP 

STKPEN 

STEPOPT 

STEPUP 


Generates the speed boundaries for the optimization search. 

Minimizes the climb and descent Hamiltonians with respect 
to speed. ' 

Controls the process of minimizing the Hamiltonian for a 
fixed energy level during climb. 

Minimizes the cruise cost for a fixed altitude and aircraft 
weight. 

Controls the process of minimizing the Hamiltonian for a 
fixed energy level during descent. 

Computes the minimum drag airspeed. 

Generates the optimum cruise cost for a given cruise weight. 

Minimizes the climb and descent Hamiltonian with respect 
to thrust. 

Calls the appropriate routine to estimate the climb fuel. 
Estima-.s climb fuel for the tri-jet. 

Estimates climb level for the twin-jet. 

Computes the value of F used as an iteration parameter to 
compute a trajectory with the desired range. 

Generates the lower (EPR) limit for climb optimization, 
and the upper limit for descent optimization. 

Prints pre-step-climb cruise performance table. 

Prints performance table, writes data on Unit 11, and 
calls for printer plots. 

Controls computation of one optimum flight profile. 

Prints cruise table summary. Calculates summed cruise time 
and distance tables. 

Controls the addition of a step climb to the optimization. 

Takes one energy step (climb or descent) using the minimized 
Hamiltonian. 

Calculates cost of combined cruise, step climb, cruise and 
descent. 

Computes a step climb in altitude and Mach number. 
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VOPTRJ Computes the fuel used, distance traveled, time taken, total 

cost and cost/n.mi. for the climb, cruise, and descent profiles. 

WATEST Estimates the landing weight, given conditions at top of climb. 

WCLST Computes final climb, cruise and descent values for fuel, 

time and distance, where the value of X is very close to 

the optimum value, 

WLEFHV Interpolates the cruise table data to relate X, cruise 

weight, energy, fuel flow rate, altitude, and ground speed. 


Under Category 2, the routines are; 


DATTRI 

DATTWN 

CDRAG 

CDRAG2 

CDRAG3 

CLIFTT 

CLIFT2 

CLIFT3 

ENGEPR 


Block data containing engine data for the tri-jet turbofan 
engine . 

Block data containing engine data for the twin-jet aircraft. 

Calls appropriate routine to compute the drag coefficient. 

Computes the drag coefficient for the tri-jet aircraft. 

Computes the drag coefficient for the twin-jet aircraft. 

Calls appropriate routine to compute the lift coefficient. 

Computes the lift coefficient as a function of Mach number, 
altitude, and angle-of-attack for the tri-jut aircraft. 

Computes the lift coefficient for the twin-jet aircraft. 

Calls appropriate routine to compute engine thrust and 
fuel flow rate. 


ENGEP2 


ENGEP3 


ENGIDL 

SPLMT 


Computes the engine thrust and fuel flow rate as functions 
of altitude, Mach number, temperature variations, and EPR 
setting for the tri-jet. 

Computes the engine thrust and fuel flow rate for the twin- 
jet aircraft. 

Com put es^tliTust^nd=^ eX-flow— rate— fox— idle— throstle-. 


Computes aircraft speed limits during climb or descent. 


TRIM Computes the thrust and angle-of-attack for maintaining 

constant speed levels for a given altitude and cruise weight. 
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Under Category 3, the routines are: 

ATLOW Generates atmospheric density, pressure, temperature, and 

speed-of-sound as functions of altitude, 

WIND Computes the wind vector and its effect along the ground 

track of the aircraft. 

WINDIN Reads in data and sets up the wind profile as a function 
of altitude. 


Under Category 4, the routines are: 

BANNER Writes the heading at the beginning of the run. 

CHEKIN Checks input quantities to be sure they are within 
reasonable ranges. 

CONDAT Block data containing values for most program constants. 

DBLSRC Performs a linear double table look-up. 

FIAS Converts indicated airspeed in feet/second to Mach number. 

FIASM Converts Mach number to indicated airspeed in knots. 

FMIN Minimizes a function by a Fibonacci search to within 1/144 

of the search interval. 

FMIN2 Minimizes a function by a Fibonacci search to within 1/34 

of the search interval. 

I CLOCK Changes time in seconds into hours, minutes, and seconds. 

JTRUNC Truncates a monotonically decreasing series from an array 
of changing values. 

LSQPOL Obtains a polynomial based on a least-squares fit to a 
set of data. 


MATTNV 
NICER 
PAGE 
PICTUR 
POLY El 

POLY 2 


Inverts a matrix. 


Generates boundaries for printer plots. 

Starts a new page of printout. 

Generates printer plots. 

Evaluates a polynomial for some fixed value of the 
independent variable. 

Evaluates a polynomial for some fixed value of the two 
independent variables. 


1 
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PRTPLT Sets up printer plots when IGRAF is greater than 1. 

PRWT Prints estimated conditions at top of climb. 

SERCHI Searches for a point in a monotonically increasing array. 

SERCHD Searches for a point in the monotonically decreasing array. 

SGLSRC Performs a linear table look-up. 

TRACIT Traces subroutine calling sequence in case of program error. 

WRITE1 Writes out the trajectory summary table. 

Three CDC subroutines are called: DATE, TIME, and STRACE. These are 

for local user convenience. 


Figures 5a through 5i show how control passes through the chain of 
subroutines. This figure, when combined with the preceeding short subroutine 
descriptions, is a short guide to the total program organization. Further 
detail may be found in Appendix B. 




Figure 5a. Top Level of Program Flow. 
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Figure 5 r _ . Program Flow for One Flight Profile from PftOFIL downwards. 
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Figure ‘it. Program Flow for Climb Optimization. 
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Figure 5g. Program Flow for Descent Optimization. 



Figure 5h. Program Flow for Near-Optimum Case. 
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APPENDIX A 

TRAJECTORY OPTIMIZATION USING THE ENERGY STATE METHOD 

The purpose of this Appendix is to summarize briefly the theoreti- 
cal background and the numerical process used in the OPTIM program 
for computing the optimum vertical profile of a turbo-jet aircraft. More 
details are given in Refs. 2-6. Reference 7 presents the principles 

upon which trajectory optimization is based. In Refs. 2 and 3, 1 

Erzberger and Lee apply these principles using the energy state 

approximation to obtain a practical, efficient means of generating 

the optimum vertical profile. OPTIM is an extension of the original 

computer code developed by Erzberger and Lee and is based on their 

methods. Its application is explained in Refs. 4, 5, and 6. 

In the following sections, the theory of trajectory optimization 
is first presented. Then, the application of this theory to minimizing 
the direct operating cost (DOC) of an aircraft traveling over a fixed 
range is outlined. This is followed by a discussion of the details of 
going from theortical expressions to a practical computer code. The 
theoretical points are presented without proof, for conciseness. The 
reader wanting more detail should review the references. 


Theoretical Principles 


In Ref. 7., a description is given of the requirements for solving 
an optimization problem involving a continuous dynamic system with no 
terminal constraints but with fixed terminal time. This description is 
repeated here because it presents the basic principles which extend to 
the aircraft profile optimization problem. 
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A system (the aircraft) is governed er js8°M^r differential 
equations 

x - f(x,u # t) ; xft^) 8 iven » (A.l) 

C o 1 C i C f 5 

where x is the n-dimens tonal state vector and u is the m-dimensionai con- 
trol vector. The cost function which is to be minimized is of the form 

J - 4(x(tf)»tf) + j* ^ L(x,u,t) dt. (A. 2) 

* o 

Here, $ is the terminal cost function# and L is the cost per unit time along 
the trajectory. The problem is to find the sequence of controls u(t) that 
minimize J. 


First, the system equations are adjoined to J with the multiplier 
vector A(t): 

J - 4>(x<t f ),t f ) + / ^(L(x,u,t)+X T (t){f<x,u,t)-x}) dt. (A.3) 


Then the Hamiltonian function is defined as 

H(x,u,t) ■ L(x,u,t) + X T (t)f (x,u,t) . 

Equation (A.3) is integrated by parts to yield 

J - ♦(*(t f ),t f ) - X T (t f ) x(t f ) + X T (t 0 ) x(t 0 ) 

+ p (HCx.u.t) + X T (t) x(t)} dt. 


(A.4) 


(A. 5) 


Next, the change in J due to variations in u(t) and x(t) is con- 


sidered for fixed t and t,j 

O L 


% 


(A. 6) 


ORir’i'v 

OF P JCis i. 


>VI 


/ 1 j(s + ;l ) 5x + s* ! dt - 


The elements of \(t) are chosen to cause the coefficients of <5x in Eq* 
(A, 6) to vanish under the integral and at t f : 

!T Z\\ Dl. \ T 3f J,. x _ 3$ 


.\ J 


Dl. V T »f x T r<K . 09 

3x “ “ " 3x • C f } " 3x * 


(A-7) 


Equations (A. 7) are called t lie co-state equations. Then, Eq . (A. hi becomes 


/ tf 

vSu dt . 
du 

t 

o 

For J to be minimum, 6J must be zero for arbitrary u(t); 
for no bounds on u, 


(A. 8) 

this implies that 


3H 

3u 


0 


t N t < t, 

o — — f 


(A.9) 


on the optimum path. If the control variables are constrained as 


C(u,t) v o, (A. 10) 

then for u(t) to be minimizing, we must have SJ ^ 0 for all admissible u(t). 
This Implies, from Eq. (A. 8) that 

|“ <5u A <5 h > 0, (A. XI) 


for all t and all admissible (Su(t). In other words, H must be minimized 
over the set of alt possible u; this is known as the minimum principle [7], 


tn summary, to solve for u(t) that minim ir.es J, the differential 
equations (A.l) and (A. 7) must be solved simultaneously, where u(t) is 
determined from F.qs. (A. 1 !) or (A. 11). The boundary conditions on the state 
x at t and \ at t ^ are specified, resulting in a two-point boundary-value 
problem. 
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(A. 12) 


Each element of Eq. (A. 12) is zero on the optimum trajectory, from which 
we can conclude that H is constant on the optimum trajectory. This latter 
point is used in the analysis presented in Refs. 2 and 3. 


Application to Aircraft Profile Optimization Using the Energy 
State Approximation 


Here we are concerned with applying the above theory to the pro- 
blem of choosing the thrust and airspeed values to control the aircraft 
vertical profile in going from one point to another. The cost function 
J is the direct operating cost (DOC) which is the sum of fuel and time 
costs. This is, in integral form, 


J 


n<C t W + C )dt 



(A. 13) 


where is the cost of fuel ($/lb), w is fuel flow rate (lb/hr), is the 
cost of time ($/hr), C d is the direct operating cost, and t f is the 
time to fly the specified distance traveled d^.. It is also assumed that the 
typical vertical profile is as shown in Fig. A.l - that is, it contains 
climb, cruise, and descent portions which have the constraint that 


where 


d + d , < d- 
up dn~ f 


(A. 14) 


up 


dn 


= x(t^) = the distance traveled from the start point to 

where the cruise segment begins (at time t - t ,). 

c i 

a d r -x(t r )= the distance traveled from the end of 
f cf 

cruise (at time t « t^ f ) to w * irre the descent segment ends. 
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Figure A • 1 . Assumed Structure of Optimum Trajoctoricc 
Thus, the cost function (Eq. (A. 13)) can be rewritten as 


(A.l 


■> -/ tCl C d dt+ W £ - d „p - V* + C d dt ’ 

0 *cf 

- C d dt + (d f - x(t cl )- [d f - x(t cf )])* + ft C d dt. 


"Cf 


illMhere ijf is the cost per unit distance while in cruise. 


Simplified point-mass equations of longitudinal motion of the air- 
craft are 


V * (T-D)/m - g sin y , 
h * V sin Y , 


<AJ 


X - Vg , 


where the flight path angle (Y) dynamics and weight loss due to fuel burn 
are neglected. Here, 


V = airspeed, (magnitude of aircraft velocity V with respect 
to the air mass), 


V = ground speed (magnitude of + V cos y), 

— « horizontal wind velocity, 

w 


T 


altitude , 
aircraft mass, 
thrust , 
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drag, and 
horizontal range. 
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Here, the range equation is based on the ground speed V (the vector sura 

© 

of the horizontal velocity of the aircraft with respect to the air mass 
and the wind velocity). Also, it is assumed that lift L ® mg cos y. 

The effect of weight loss is accounted for by continuously updating weight 
without adding another state variable. 

The objective of this development is to simplify the optimization problem so 
that it can be solved in an on-board computer. This is done by use of the 
energy state approximation which is now presented £8]. Specific energy E is 
defined as 

E = h + V 2 /2g, (A. 17) 

which is the sum of potential and kinetic energy per unit mass. Its 
time derivative is found to be 

E « V (T-D) /mg . (A. 18) 


The energy state approximation is based on the assumption that potential 
and kinetic energy can be interchanged instantaneously. In this approxima- 
tion, the energy state variable replaces altitude and airspeed state 
variables 18). Thus, Eq. (A. 17) can be used in place of V and h in Eq. 

(A. 16). 


It is assumed that the aircraft specific energy increases monotonicai ly 
during climb and decreases monotonically during descent. This assumption 
is used in the development to change the independent variable in Eq. (A. 15) 
from time to energy. This uses the transformation 


dt 



(A. 19) 


It is mathematically convenient to evaluate the last integral in Eq* 
(A. 15) backwards in time so that the energy state is monotonically increas- 
ing during its evaluation. This means that the running distance (range) 
variable during the descent can be measured backwards from the end point. 
Thus we can think of range measured in two ways as shown in Fig. A. 2. 
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In this sketch, 

) * range measured on the way up in forward time t, 
j) 5 * value of when initial cruise is reached, 

sured on the way up in backward time t , 
x dn^ T cf^ = value x dn when final cruise is reached ( T c f“ I t c f~ t f I ) ■ 

Also, we define the variable x to be range traveled during climb and 
descent. The distance traveled during cruise is then constrained to foe 
(df ~ x ) • We can then see that an incremental change dr in the range 

dn* 


W 0 “ 

range 

up (t ci>“ 

value 

CL 

H 

II 

range 

d„ <T cf>- 

value 


variable x is equivalent to incremental changes in both x and x 

up 

That is, 


dx 


d(x + X . ) a 
up dn 


(A. 20) 


t rom this discussion, the second of Eqs. (A. 15) can be written as 

J - P cic d dt + (d f - x up (t tl ) - x dn (T cf ))^ + r Cf I C d I dx . 

n J 


We use Eq. (A. 19) and the transformation 


di 


dE 

m 


(A. 21) 


(A. 22) 
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to rewrite Eq. (A. 15) as 


■m «) + « £ -<v E d> + x d» <E of ,)) * 

E i fe>0 

•/■•'fe-L 


(A. 23) 


Here, E^, E ci , E f , and are the values of specific energy evaluated 
at time t equal to 0 and t ci and time t evaluated at t c £ and respec- 
tively. 


Note from Eq. (A. 23) that the range x only appears as the sum of 

climb and descent distances (x + x. ). Thus, the state equation for 

up an 

this system of equations can be written as 


. d(x + x . ) 

dx _ up dn* 

dE 


dE 


E>0 1 1 E<0 


(A. 24) 


Here, V and V ^ are the equivalent ground speeds for climb and 
descent. Then, analogous to Eq. (A. 4), the Hamiltonian is 



/c.\ 

/C. \ ( 

/V \ 

/V „ \ 1 

i 

H = 

~ + 

( m ) + ^ > 


l-sis.) 1 

I 


; e 4 >o 


A E ' E >0 

'i*i Li 

( 


(A. 25) 


This can be divided as 


H = 


C , + A (V ) 
A _£HE_ 


C d + » <'W 

1 

tx)-» 

1 

1 

1 

PI- 

V 

o 

H 

|fc| 


(A. 26) 


E<0 


Now, analogous to Eq. (A. 7), the costate equation for A can be written as 
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•'K 


d Y 


3.11 

d ( x + x t ) 
up iln 


(A. 27) 


and t com Kqu . (A. 7) and (A. 21), this cost ate has tin* Unai value 


u, , . wp , . .'"fVV. 1 .., 

u '- |> '-'f a <V»„,,' »‘VW 


(A. 28) 


where is tlie erulse cost per unit distance. 


Mote, this problem could be placed in a slightly more conventional term 
by divULtnn It Into two problems - one tor climb and one-half ol the 
cruise distance and the other for decent and the other half of the cruise 
distance. Then Kqs. (A. 27) and (A. 28) would be replaced by 


OA 3H n 

» ■ ■ 


'<V 


8 «V 2 - V 1 


t*X 


up 




'».2«> 


for climb. For descend. 


;n 

3K 



0 


»(E cf ) 


jadj/2 - » to i *<E C( » 


3, 


dn 


- - * t* cl ) 


(A. 30) 


This allows splitting the Hamiltonian defined In Kq. (A.2fc) and allows for 
X(E t>l ) i X(E f ). In fact, in the* actual Implement at ion E f t i K cf because 
optimum cruise energy changes as lucl is burned off. The principal results 
are unchanged, however. — — ^ 

11ms, from Kq. (A.1U, <A.2‘») and (A. 10) the trajectory opt Ira l sat Ion 
problem becomes 
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(A. 31) 


Thus, the optimization problem reduces to solving pointwise minimum 
values of the algebraic functions defined by Eq. (A.31) during the climb 
and descent portions of the trajectory. 

Equations (A. 29) and (A. 30) are the transversality condition for the 

free final state problem (d + d , < d c ) with terminal cost. Thus, the 

up an t 

constant value c' A~“for climb and descent is found to be the negative of 
the cost per unit distance for cruise. 

The cruise cost (« - X) is found by assuming that the aircraft is 
in static equilibrium during cruise (T = D) , and that 

? ft) • «&£ 

V + V I . 

C W 1 

In other v rds, for any cruise altitude, there is an optimum thrust and 
airspeed V c such that the cost per unit distance ^(E c ) is minimized. The 
optimum cruise cost as a function of cruise energy is typically of the 
shape shown in Fig. A. 3. Thus, there is also an optimum cruise energy 
E copt w ^ ere crui - sc cost ^(E copt ) is minimized. If the range is long 
enough so that there is sufficient range to reach optimum cruise energy 

E copt durin § climb, it should be done, and the cruise conditions should be 

set so that x(E ) =-tL(E J. 

c r ' copt 


where 


*<= c > 


gc 
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Figure A. 3. Optimum Cruise Cost as a Function of Cruise Energy 

For the case where there is no cruise segment (d* * d + d, ). 

f up dn' 9 

the cost function contains only integral terms* Then, the transversality 
condition yields X- - ^(t c ). That is, \ would be the negative of 

where Mt c ) is the optimum cost for cruising at the highest point 
reached on the climb trajectory. 

The optimum cruise energy % CQpt is only specifically reached when 
there is range enough to climb to and descend from the optimum altitude/ 
airspeed values, where (E c ) is minimum. For ranges less than this 
value, the maximum value of E c that is reached is a free variable less than 
the optimum value. Its choice is made to optimize the cost function of 
Eq. (A. 23). 


From Eqs. (A. 23) and (A. 25), one can write 

M. ... T af(d f - d u p -“W *< E e»l . 

3E " H + Tf ■ 0. 


at E * E . This is 
c 
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where d is the cruise distance* and H Is the total value of H (H + H . ) 
c c 'up dn' 

at the cruise point. Thus, Eq. (A. 14) can be used, along with other 
characteristics ct ip and H, to determine the relationship between E » and 
d c# The Hamiltonian evaluated at E * is the cost penalty to achieve a 
unit increase in cruise energy. For N 0, Eq. (A. 34) can be written as 

- -H c /(3*/3K) k . F (A. 35) 

V 

Figure A. 4 shows the family of trajectories which have this characteristic. 

Those occur at values of R below E where Htp / n K * 0 (see Fig. A, 3). 

c copt 

That is, non-zero cruise segments occur at short ranges with cruise energies 
less than the optimum energy value for long range. 


For the case where H • 0, Is zero for 3^/3E < 0. The distance 

c c 

d c can be non-zero only at optimum cruise energy where 3^/3E ■ 0. This 
family of trajectories is shown in Fig. A. 5. 


Thus, we have a situation where positive values of H dictate one 

c 

type of trajectory and zero values dictate another. In Ref. 2, it is 
shown that if the aircraft engine specific fuel consumption S_„ is in- 
dependent of the thrust T (so that w • Sp^T) , then the structure of the 
trajectories will be like Fig. A. 5 with no cruise segment occuring except 
at E CO p t * (This implies that the Hamiltonian H c is zero at the maximum energy 

point). For this case, the optimum thrust setting for climb is T , and 

max 

the optimum setting for descent is T, ,, . 

idle 

If the engine specific fuel consumption is dependent on thrust, and 
the thrust values are not constrained during climb or descent, it is 
shown in Hof. 2, that the Hamiltonian H is again zero at the cruise 
energy, and again the trajectory structure is like those of Fig. A. 5. 

18 dependent on thrust for the aircraft models used in the OPTIM 
program. 
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If the S pc is dependent on thrust, and constrained to the maximum 
value for climb and to the minimum idle value for descent, then the 
Hamiltonian is positive at cruise. This causes positive cruise segments 
according to Eq. (A. 35) at cruise energies below the optimum. For this 
case, the optimum trajectories will have shapes similar to Fig. A. 4. 
These trajectories are slightly less efficient than those of Fig. A. 5. 
because one less control is available for optimization. 

Some Mechanization Details of the Computer Program 


The remaining sections of this Appendix describe how the previous 
theoretical material has been utilized to construct an offline computer 
program for generating optimum vertical profiles for models of two turbo- 
jet aircraft provided by NASA. This material is presented in an alternate 
way in Ref. 3, and the program is referred to here as OFTIM, 

By examining the specific fuel consumption data of the turbojet engine 
models, it was determined that S FC is dependent on thrust. Thus, two 
types of short range profiles must be considered - those represented 
by Fig. A. 4 (Type 1 profile) when thrust is constrained and airspeed is 
the single control - and those represented by Fig. A. 5 (Type 2 profile) 
when both thrust and airspeed are used as controls. 

The solution to optimum climb and descent profiles is found by 

minimizing the Hamiltonian expressed in Eqs. (A. 31). The independent 

variable (energy) Is stepped along In fixed increments (e.g., 500 ft), 

and the Hamiltonian is minimized at each energy setting. Minimization 

occurs by finding the best values of airspeed (V , V^) an< * P° &s ibly 

thrust (it , tt, ) so that the climb function and the descent function are 
up an 

Individually minimized. 

To solve Eqs. (A. 31) requires knowing two more quantities: 

\KE ) - the cruise cost per twi* distance. This comes from evaluating 
c Eq. (A. 32) at the deal ced cruise altitude. 
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E - the cruise energy. This is a function of the cruise altitude 
and the associated cruise airspeed obtained in Eq. (A. 32). 

Note that for the Type 2 profile at short ranges, there is no cruise 
segment. In this case, the maximum energy achieved at maximum altitude 
is referred to as the cruise energy E^. At that altitude, there still 
is defined a minimum cruise cost according to Eq. (A. 32). 

For the Type 1 trajectory of short range, there exists a non-zero 
cruise segment which is determined by use of Eq . (A. 35). To solve Eq. 

(A. 33) requires that the Hamiltonian defined by Eqs. (A. 31) be solved 
at the point of transition from climb-to-cruise. It also requires 
knowing the slope ?ii//8E of the cruise cost for a change in cruise energy at 
that point. 

Cruise Optimization 


The first step that must be taken to compute optimum trajectories 
is to derive the optimum cruise cost ^ and its derivative 3^/3E. This 
is done by computing what is referred to as the "cruise table* 1 . The 
parameters that affect this table are the assumed cruise weight, the 
wind profile, and the lift L, drag D, thrust T, and fuel flow rate w 
characteristics of the aircraft. The optimization process searches over 
the acceptable ranges of altitude and airspeed for a given weight. The 
results are collected in tabular form for a series of different assumed 
cruise weights. 


't . 

5r 

fe 



Again, the minimum cost of flight during cruise per unit distance 
for a fixed cruise weight is found by 




min 

V 


c f w + G t 


gc 1 


(A. 36) 
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This assumes that the aircraft is in static equilibrium during cruise, i.e., 

T cos a * D f (A. 37) 

L + T sin a « W, 

where the angle-of-attack ot is found by solving these equations simul- 
taneously* The altitude is stepped in 1000 ft increments from 10,000 ft 
to ceiling altitude (where maximum thrust just balances drag). At 
altitudes below ceiling altitude, the airspeed - dependent drag curve 
crosses the maximum thrust curve at two points (V^ and as illustrated 
in Fig. A. 6. Thus, for' each altitude level, the values of and 
are determined, and then ^( W C » E C ) is minimized with respect to airspeed 
V between these two limits. Restrictions are that be greater than 0.1 
Mach and that be less than 0.89 Mach or 0.84 Mach for buffet constraint reasons. 

After the cruise cost is minimized at each discrete altitude level, 
these numbers are stored in a table with altitude as the independent 
variable. Typical results are plotted in Fig. A. 7. Presented here are 
also the optimum cruise Mach number M opt and the optimum thrust setting 
EPR Qpt . After results are obtained in steps of 1000 ft, the minimum 
cost point is found as a function of altitude. In the 0PTIM program, 
the cruise table optimization results are obtained by using a Fibonacci 
search with eight Fibonacci numbers. 

The cruise table results are obtained for cruise weights varying as 
dictated by the program input. Usually, the cruise weight is incremented 
in steps of 3000 lb. Up to ten values of cruise weight can be used. For 
each cruise weight, the optimal cruise altitude, cost, speed, power 
setting, fuel flow rate, and specific energy are computed. An example 
of optimum cruise cost as a function of cruise weight is shown in Fig. A. 8. 


C limb Optimization % 

After the cruise tables are generated, the program precedes with 
obtaining tho optimum climb trajectory. This requires guessing what the 
cruise weight will be, based on the takeoff weight. This guess is 
used to obtain a trial value for (or X ) in the Hamiltonian from 
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Figure A. 8. Optimum Cruise Cost and Cruise Altitude as Functions of Cruise 
Weight for the Tri-jet Aircraft Flying into a Particular Head Wind. 


the cruise tables. The procedure to obtain this guess is based on an 
empirical formula which iterates until convergence is made. 

The climb optimization process starts by assuming X (E c ) * 1,01 ^ c ^ E C opt^ 

or 1.0 ti/ (E ) where (E ) is first obtained by setting the initial 
T c copt * c' opt 

cruise weight W ^ equal to the takeoff weight (an input) . The appropriate 
cruise tables are used to interpolate to find the corresponding value of 
E associated with 1.01 ip or 1.0 ip . Then, empirical equations are used 

C C c 

to obtain an approximation to the fuel burned to reach E c * For example, 
for the tri-jet model, the form is 

K up = °- U (E ci - E i> (1 + A - 7 W W ci /W ref < A - 38 > 

Here, is the takeoff aircraft energy, ^ re f is a reference weight (136000 lb) 
for the tri-jet, and W ci is the previous value of cruise weight. Then, 
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Che cruise weight is updated at W 


ci 


W . - F . 
ci up 


This process is 


repeated until the difference in consecutive estimates of F falls below 

up 

100 lb. 


When the cruise weight estimate is obtained, the corresponding values 

of E and X(E ) are obtained from the cruise tables. Then, the program 
c c 

is ready to generate points on the optimum climb trajectory. This is done 
by stepping along at discrete increments of specific energy and minimizing 
the Hamiltonian function 


H 

up 


(E) 


c f w + C t - X(E C ) (V + V w ) 
E 


(A. 39) 


at each point. (This is the first of Eqs. (A. 31)). That is, the program 
starts with initial energy E q * h Q + V q /2g. It steps the energy a 
fixed amount AE (say 500 ft). At this point, it searches over airspeed 
V (and possibly thrust setting ) so that Eq. (A. 39) is minimized. For 
the turbojet engines, thrust is governed by EPR settings which vary 
between 1.1 (idle thrust) and some maximum value less than 2.4. The air- 
speed has an upper limit governed by 

a) . 0.89 or 0.84 Mach (buffet limits) , 

b) . 250 kt (IAS) below 10,000 ft for ATC restrictions, 

c) . / 2g(E-h) which insures that the aircraft climbs, and 

d) , the upper value shown in Fig. A. 6 where max thrust 

equals drag. 


The lower limit is governed by 

a) . , the lower value shown in Fig. A. 6 where max thrust 

equals drag, __ _____ 

b) . 0.1 Mach 

c) . 5 ft/sec less than the previous value of V to limit large 

jumps in flight path angle. 

The Fibonacci search technique is again used to determine V and tt 
which minimize Eq. (A. 39) for the fixed value of energy E. The value 
chosen for aitapecd Is accura t€L_tjQ_wlthin — 0056 Mach, and EPR is accurate 
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energy rate E (Eq. (A. 16)) and altitude h: 


h = E -V 2 /2g 
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(A. 40) 


From these, approximate values of time, range, flight path angle, and 
fuel burned are obtained from 


At = AE/E, (A. 41) 

sin y * (Ah/At)/V, 

x * ^Ax ; Ax = At 

F = >:AF ; Af = w At. 


The ground speed is found as 

V g = V w cos(* g - * w ) + y V w 2 cos 2 W g - * w ) + V 2 - V w 2 . (A. 42) 

where ij; and ^ are the desired aircraft ground heading and wind heading, 
g w 

respectively. 


The above process is repeated by stepping along energy in incre- 
ments of Ae until E c is reached. The last value of Eq. (A. 39) is 
stored for possible use in evaluating the cruise distance. 

The above climb optimization procedure is repeated with x set to various 
multiples of the optimum cruise value ^ until the total range of flight con- 
verges to the appropriate value. This is discussed in further detail later. 


Descent Optimization 

The descent optimization is very similar to the climb optimization 
with regard to the equations which are evaluated. The optimization pro- 
cess requires estimated values of and at the beginning of descent, 
and an estimate of weight at the end of descent. The method used 
to obtain these estimates is discussed in the next section. 

If there is a cruise portion of flight, fuel will be burned off during 
cruise. Thus, the value of E^, 4‘ c > and W c ^ at the beginning of descent 
will be different than at the beginning of cruise. If there is no cruise 
portion, then these values will be identical. 
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The descent profile is obtained by starting at the final energy 
state and then going backwards in time. The energy rate is constrained 
to be negative with respect to forward time. 


Similar descent profile constraints exist on airspeed as for those 
of the climb profile. The thrust level is on or near the idle value 
during descent. 


Cruise Fuel Burn 


To estimate the final weight during cruise (W^) and landing (W^ 
the following steps are taken: 


1) . Determine ^ , the initial cruise cost based on the initial 

c 

cruisjc_we±ght obtained from the climb optimization. 

2) Use the initial cruise weight to compute the fuel flow 
rate w(^ c ) from the cruise table. 

3) . Estimate the cruise range by use of empirical equations; i.e., 

for the tri-jet model: 

P = ^c^copt * 1#0 or 1 - 01 > (A.43) 

d = b.P 4 + b,P 3 + b,P 2 + b,P + b, . 
cl 2 3 4 5 

4) . Compute the cruise f^el as 

F c - V v (A,44) 


5). Estimate the average cruise weight as 


W « W ~ 0.5F 
c c c 


6). Use the cruise table to obtain the corresponding cruise cost 
• 

^ c t altitude h c> fuel flow rate wtfp ), airspeed V f and wind 


speed v w (h) . 


7). Recompute Eq. (A. 44} and then find the final cruise weight. 


W = W - F . 
cf ci c 


8). Use the value W .. in the cruise tables to obtain ^(W £ ). 

c r ct 

As with the climb, set X = 1.0 i|» c (W cf ) or 1.01 ^ c (W cf ). 


67 


ORIGINAL PAGE IS 
OF POOR QUALITY 


9). Use this value of > to ootain h cf and E( . f from the cruise 

tables. 11, ese are the end conditions for the descent trajectory 
obtained backwards in time. 


10 ). 


hstlmato the landing weight from empirical formulas^ e.g., 
the tri-jet: 

P - 1.01 or 1.0 


w r = w rf - 


+ C 2 P 




for 

(A. 47) 


The values ot X, K^j. and Wj. obtained by the above procedure are used for 
obtaining the optimum descent trajectory. The descent portion of the 
Hamiltonian is of the form 


w + 


dn 


(10 = 


C t - X(K el> 


<V 


(A. 48) 


this f urn * ion 


is also minimi zod at on ch of 


t ho ivt'n values of onor>;v. 


Aftei the tirst descent profile is completed, a new estimate of 
cruise distance is obtained by using Kq. (A. 35), or 


d 

c 



+ H dn )/(at|»/3K) 


(A. 49) 


Tiien, step (4) above is repeated to obta-in— an— improved cruise fuel burn. 
Then, the improved landing weight estimate is 


W f = W i 



F 

c 


+ K dn ) 


(A. 50) 


Hie landing trajectory is reoptimized with this now value of 
landing weight. Then, improved values of total range traveled, time 
required, and fuel burned during climb, cruise, and descent are made. — 


For short range flight, the above steps assumed that a Type 1 trajec- 
tory is generated because thrust is const ra+ned-te-max imurn value during 
climb ar.d idle value during descent. If thrust is free, then a Type 2 

result . with no mils'* pot t ion. For this ease, the steps 
required to estimate cruise distance ,1^ and final cruise cost, weight, and 
energy can be eliminated. 
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Cruise Cost Estimation 


The first climb and descent profiles are generated with X « ^c^ E copt^ 

for free thrust or X « 1.01 ^ C ( E C t ) £° r fixed thrust. The resulting 

range of the optimum profile has a given value defined as R . An example 

is shown for X * 1*01 (E . ) in Fig. A. 9* If the total range desired 

c opt 

is greater than R , then the climb associated with X * (E ) is 
max T c copt 

used, the required cruise distance is computed, and the descent profile 
is recomputed to produce the correct overall range. 


If the desired range is less than R , X is next set to a value 

max 

corresponding to a cruise altitude just over 10,000 ft. (This typically 
has the value X ■ 1.3 - 1.5 <P c ) . The optimum profile is recomputed, and 
the associated range is referred to as R min » (See Fig. A. 9) 

If the desired range is between R . and R , then an iterative 

min max 

process is used to obtain i|>( E c ) and the associated desired range. Itera- 
tions are stopped when the total range traveled is within some e of the 
desired range. (In OPTIM, e is set at 5 n.mi.) 

If the cruise altitude is fixed, or a step climb is used between two 
fixed altitudes, then X is set to ^( E copt ) corresponding to the fixed 
altitude^ In this ca e, no iteration on X is required, 
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APPENDIX B 


OPTIM SUBROUTINE DESCRIPTION 


This section contains an explanation of the data storage for program 
OPTIM. Following that is an explanation of the main program, its two 
principal subprograms, and then the remaining subroutines and functions 
in alphabetical order. 

Data Storage 

The major part of the data communications between subroutines in 
OPTIM is through labelled common statements. There are ten such commons. 
Their names and a short description of each are: 

CCDE Cruise, climb, descent variables. 

CONST Constants. 

CRUISE Cruise table and associated variables. 

DESCRP Assorted variables. 

ERROR Error traceback information. 

GRAPH Data to be written to Unit 11 and associated variables. 
Includes the final climb and descent trajectories. 

INPUT All input parameters* 

TOA Time-o f-arrival and step climb variables. 

TRIJET Engine data, tri-jet. 

TWINJT Engine data, twin-jet. 

WINDP Wind input data and associated parameters. 

As a convenience * the CDC UPDATE capability is used to insert COMMON 
statements into source decks. This facilitates changing items in COMMON 
v ith no loss of program portability, because UPDATE produces a compile file 
which is directly iistable, editable, and compilable by any standard 
FORTRAN. 

Subroutines 

Following are more detailed descriptions of the main program OPTIM 
and its subroutines and functions. 
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MAIN PROGRAM: OPTIM 



This program synthesizes a fixed range, minimum fuel or direct operating 
cost trajectory. The overall process is explained in Section IV and 
Appendix A. The main program calls a subroutine (ALLIN) which reads all 
required input, and then calls the principal control subroutine (0PTM56). 

Subroutines called: Commons used: 

ALLIN None 

OPTM56 


OPTM56 

QPTM56 controls the main program computations. It first calls 
CHEKIN to check that all input quantities are within reasonable limits. 

It may then call OPTTOA to execute a run with fixed time-of- arrival, if 
ICTAB * 2. Otherwise, it calls in sequence PROFIL to compute the optimum 
requested profile, STEP to adjust the trajectory to include a step climb, 
if appropriate, and PRFTBL to output the results of the optimization. If 
any error has ocurred during the computation, TRACIT is called to output 
the traceback information. 


Subroutines called: 


Commons used: 


CHEKIN 

OPTTOA 

PRFTBL 

PROFIL 

STEP 

TRACIT 


CCDE 

CRUISE 

DESCRP 

ERROR 

GRAPH 

INPUT 


f 
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OPTTOA 


This subroutine serves as the program executive when synthesizing a 
vertical flight path which requires minimum fuel to achieve a fixed range 
in a fixed time-of-arrival. The logic is similar to that of OPTM56 
except the following items are changed: 

1. An outer loop structure is mechanized to iterate on TC to 
achieve the desired time-of~arrival TEND. This is illustrated 
by the block diagram in Fig. 4, Chapter IV. TC is initially 
set to zero, and the minimum fuel trajectory is generated. 

The time this trajectory takes, TTIME^, is compared to TEND, 

and the subsequent logic depends on whether TTXME is greater 
or less than TEND. 

2. If TTIME q is less than TEND, a special trajectory is next gene- 
rated with the cruise speed set at that value where minimum fuel 
rate (max/L/D) is achieved. This represents the upper bound on 
length of flight time that is practical without using path 
stretching or a holding pattern. If this trajectory takes TTIMEq, 

then TEND is compared to TTIMEq. If TEND is greater than TTIHEq, 

the program stops. If TTIME o < TEND < TTIMEq, the cost of time 

TC is negative, and iterations on TC continue until TTIME is 
within 10 sec of TEND, 

3. If TTIME q is greater than TEND, the cost-of-time TC is positive. 
Iterations on TC continue until TTIME is within 10 sec of TEND, 


The logic uses the previous three values of TC to fit a quadratic 
curve to TTIME as a function of TC. The desired value TEND is then used 
to obtain the next trail value of TC. The program stops if convergence 
has not been achieved after nine trials. 


Subroutines called: 


Commons used: 


MATINV 

PRFTBL 

PROFIL 

STEP 


CODE 

CONST 

CRUISE 

DESCRP 

ERROR 

INPUT 

TOA 

WIN DP 


ALLIN 


ORIGINAL PAGE IS 
OF POOR QUALITY 


ALLIN reads all input, including the cruise table if necessary. It 
calls WINDIN to read the wind in. ALLIN also outputs the input data and 
lists the options chosen. 


Subroutines called: Comons used: 


BANNER 

CCDE 

PAGE 

CRUISE 

SERCH1 

CONST 

TRACI T 

DESCRP 

WINDIN 

ERROR 


INPUT 


WINDP 


ATLOW 

2 4 

This subroutine generates the atmospheric density (in lb sec /ft ), 

2 

atmospheric pressure (in lb/ft ), atmospheric temperature (in degrees 
Kalvin) and speed of sound (in ft/sec) at a given altitude below 20,000 
meters (65,617 feet). It also makes the appropriate modifications in 
atmospheric density and speed-of-sound to account for variations in 
standard day temperature (represented by the input DTEMP) . The 1962 
standard atmosphere is used. This version of the program does not 
calculate a new atmosphere when called &t successive times at the same 
altitude. 


BANNER 


original 

Of POOR 

BANNER is a subroutine used to write the title, date, and time of day 
at the beginning of a run. 


pace 19 
Quality 


Subroutines called: 

TIME } CDC-supplied 


Common used: 
None 


BLOCK DATA - DATTRI 

This block data contains the engine data used with the tri-jet aircraft 
model. Three tables are used to describe idle thrust, idle fuel flow, and 
maximum continuous engine pressure ratio (EPR) . 

Subroutines called: Common used: 

None TRIJET 


BLOCK DATA - DATTW N 

This block data contains numerical characteristics of the turbofan 
engine used with the twin-jet aircraft model. Seven tables are used to 
describe idle thrust and fuel flow for bleed valves open and closed, 
altitude of surge bleed valve closure, maximum EPR for climb and cruise, 
and Mach number corrections. 

Common used: 

TWINJT 


Subroutines called: 


CDRAG 


This subroutine calls the appropriate routine to compute the aircraft 
drag coefficient CD based on the particular aircraft model selected by 
the input variable IAC. Currently, two models are available, but logic 
is present to use up to four different aircraft. 

Subroutines called: Commons used: 

CPRAG1* None 

CDRAG2 

CD RAG 3 

CDRAG4* 

* not included with program 


CDRAG2 

This subroutine computes the drag coefficient CD for some given 
Mach number EM and lift coefficient CL for a medium range tri-jet trans- 
port aircraft model. The value is computed from the coefficients of a 
polynomial stored in the array COEFF, 

Subroutines called: Commons used: 

POLY 2 None 


CDRAG3 

This subroutine computes the drag coefficient CD for some given 
Mach number EM and lift coefficient CL for a medium range twin jet trans- 
port aircraft model. CD is computed by polynomial evaluation, including 
interpolation of the polynomial and its first derivative in certain 
regions, as necessary. 

Subroutines called: Commons used: 

POLY 2 None 
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CHEKIN 


CHEKIN tests several of the input quantities* such as the weight 
and time inputs* to ensure that they are within reasonable limits. 

Subroutines called: Commons used: 

None INPUT 


CLIFTT 

CLIFTT calls the appropriate routine to compute the aircraft lift 
coefficient for the particular aircraft model selected by the input 
variable IAC. Currently* two models are available* but logic is present 
to use up to four different aircraft. 


i . 

ft .. 

I 

& 


Subroutines called: 

CLIFT1* 

CLIFT2 
CLIFT 3 
CLIFT4* 


^Commons used: 


None 


* not included with program. 


CL I FT 2 


This subroutine computes the lift coefficient CL for a medium range 
tri-jet transport aircraft as a function of Mach number EM, altitude H, 
and angle-of-attack ALPHAP. The lift coefficient consists of three terms: 

C L " °L (ba8ic) + ^0 + C La “ 

The first term (basic) is a polynomial function of angle-of-attack a. 

The value of this term is checked against the buffet boundary expressed 
as a polynomial of Mach number* The second terra is a polynomial of 
Mach number with altitude as the parameter. The third term is also 
a polynomial of Mach number. The coefficients of the polynomial are fit 
for different altitudes. 

Subroutines called; Commons used: 

POLY t ■■ None 


CLIFT3 

This subroutine computes the lift coefficient CL for a medium range 
twin jet transport aircraft as a function of Mach number EM, altitude H, 
and angle-of-attack ALPHA, The lift coefficient consists of three terms: 

Cl = C L (baslC> + C Lo +C La° 

The first term (basic) is a function of angle-of-attack. The second 
term C_ is a function of altitude and Mach number. The third term C T 
is also a function of altitude and Mach number. These terras are determined 
by table lookup. 

Subroutines called: Commons used: 

DBLSR C None 

•~S^RCtt£== — 
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CLIMB 


This subroutine controls the computation of the climb portion of 
the fligh*. profile. CLIMB calls FULEST to estimate the fuel use during 
climb. It then uses the associated X from the cruise table computation 
to optimize the climb trajectory. The optimization is calculated in 
FCLMB6 at a series of energy steps. The energy step size is computed 
in ESTEP, and the incremental time, distance, fuel used, and associated 
quantities are computed in STEPEN at each step. All climb variables 
are stored at each step for later output and, if desired, plotting. 


After the climb is completed, the actual fuel use is found. If 
the estimated and actual fuel use differ by more then 200 pounds, the 
climb is recomputed with the X associated with the new weight. Only 
one iteration of climb is permitted. 


Subroutines called: Commons used: 


ATLOW 

CCDE 

ESTEP 

CONST 

FCLMB6 

DESCRP 

FIAS 

ERROR 

FULEST 

GRAPH 

I CLOCK 

INPUT 

PRWT 

WINDP 

SGLSRC 


STEPEN 


VILEFHV 



CRUISR 

CRUISR is used in the step climb option only. It is assumed that the 
aircraft cruises from weight to W 2 (with being greater than W 2 , in 
the initial pre-step climb cruise table). CRUISR returns cost, energy, 
fuel use rate and Mach number at the new cruine weight It also 

computes distance travelled in going from to W 2 * 

Subroutines called: Commons used; 

SERCHD CCDE 

CRUISE 

DESCRP 

ERROR 

TOA — 
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CRUISX 


CRUISX looks up in the cruise table, and returns the flight para- 
meters associated with crusing for a horizontal distance of RANGEX, 
starting at weight W y On later entries, it computes the effect of 
crusing an increment of RANGEX from the previously computed position 
(without looking up a new weight). 


Subroutines called: 


Commons used: 


CTABLE 

SERCHD 

SERCHI 


CODE 

CRUISE 

DESCRP 

ERROR 




CRZ0P5 


ORIGINAL PAGE IS 
OF POOR QUALITY 


This subroutine generates the cruise table. The cruise table includes 
the minimum cruise cost if/(W) » altitude h, and airspeed V as functions of 
weight W. For each given weight W, the optimum cruise cost, altitude and 
speed are computed as follows: 


For some given altitude h: 


a) determine the minimum drag speed V^, 


b) compute the maximum thrust T available for this 
speed and altitude, m 


c) if T is greater than D and V proceed; otherwise 
' max mD 


the ceiling has already been reached, and proceed to 2. 

d) compute the lower and upper permissible speeds and 


V., where the maximum thrust curve intersects the drag curve, 

D 


e) minimize the following cost function, 


C f £ + C t 


V V b 


where 


the fuel cost in $/lb, 
fuel flow rate in Ib/hr, 
time cost in $/hr, 
ground speed in ft/sec. 


Steps (a)-(c) are represted for all permissible altitudes 
ranging from 10,000 ft to either HMAX or the maximum ceiling. 

For a fixed altitude run, values are calculated for only the 
input altitude. For fixed TOA, the minimum altitude is 20,000 ft. 


2. Determine the minimal cruise cost X* (W) and its associated h 


and V as follows: 


A* (W) = min \p 

h 


(For a fixed altitude run, this minimization step is omitted.) 


Steps (1) and (2) are repeated for different weight values starting with 
WEIGHT and decreasing in steps of DEW through WN. 
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During a fixed time-of-arrival run, on the second iteration (IMFD-1) , 
this subroutine performs the same process to determine the airspeed at each 
altitude and cruise weight where fuel flow rate f is minimized. This 
occurs at near maximum lift/drag. 


Sri- 



The output of CRZ0P5 may be written on Unit 8 which can then be 
stored for later use. The output is also printed cn Unit 6. 



j 


ja-y 


g:;1 

Subroucines called: 

Commons used 


FIASM 

CCDE 


FMIN2 

CRUISE 

BJr 

FOPT 

CONST 

i?y- 

JTRUNC 

DESCRP 


LSQPOL 

ERROR 

- \i ■ 

PAGE 

INPUT 


SERCHI 

VJINDP 



TOA 


CTABLE — 

CTABLE looks up and returns the set of parameters in the cruise 
table corresponding to weight WCRUZ. 


Subroutines called : 
SERCHD 


Commons used: 

CODE 

CRUISE 

DESCRP 

ERROR 



% 
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DBLSRC 


ORIGINAL PAGE fS 
OF. POOR QUALITY 


This function performs a double table lookup. Given a function f(x,y), 
this function interpolates the appropriate arrays to obtain approximate 
values of f (A f B) . The four points which surround (A,B) are first found* 
and the function is evaluated at each. Then these values are interpolated, 
first on x and then on y, to obtain the approximate solution. 

Subroutines called: Commons used: 

SERCHI None 


DESCND 

DESCND controls the computation of the descent portion of the flight 
profile. It first calls WATEST to estimate the remaining cruise rang e , 
cruise fuel, and descent fuel. After this point, the computation is 
similar to that of CLIMB. ESTEP specifies energy step increments, at 
each of which the trajectory is optimized within FDESCN6 . STEPEN is 
called to calculate additional dependent parameters at each step, and 
all quantities are stored for later output and retrieval. In addition, 
DESCND computes cumulative descent cost and total cost (climb plus 
descent) at each step. 

If DESCND has been called during a constrained descent (DESPC) , it 
saves in COMMON the value of the altitude and fuel use at which the 
unconstrained descent crosses the critical altitude, HCABSL . 


Subroutines called: 


Commons used; 


CTABLE 

ESTEP 

FDSCN6 

I CLOCK 

PRWT 

SERCHI 

STEPEN 

WATEST 


CCDE 

CONST 

DESCRP 

ERROR 

GRAPH 

^nput 

WINDP 
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DESPC 


This subroutine modifies the final portion of the cruise and descent 
profile to include the sink rate constraint HDOTC from cruise altitude 
down to altitude HCABSL. DESPC first determines where the descent crosses 
HCABSL. It then generates a new, constrained upper portion of the descent 
with a ramped Mach number and constant sink rate; it modifies the descent 
table accordingly. Finally, DESPC calls VOPTRJ to join climb and descent 
with a new cruise segment. 


Subroutines called: 


Commons u^ed: 


ATLOW 

DESCND 

FIASM 

FOPT 

I CLOCK 

TRIM 

VOP 

VOPTRJ 

WIND 


CCDE 

CONST 

CRUISE 

D^SCRP 

GRVPH 

ERROR 

INPUT 

WINDP 


DRAGC 

This function computes the drag coefficient as a function of 

Mach number and lift coefficient C^, using subroutine CDRAG and the_weight, 

aircraft and atmospheric parameters previously stored in common. It 

assumes lift equals weight to compute C ♦ 

Li 

Subroutines called: COMMONS Used: 

CDRAG CCDE 

DESCRP 

INPUT 
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ENGEPR 


This subroutine calls the appropriate routine to compute the aircraft 
maximum thrust and EPR, the thrust associated with the input EPR, and the 
fuel flow rate* The engine model is associated with the particular air- 
craft model selected by the input variable IAC. Currently, two models 
are available, but logic is presented to use up to four different aircraft 

Subroutines called: Commons used: 

ENGEP1* None 

ENGEP2 

ENGEP3 

ENGEP4* 

* 


not included in program 



ENGEP2 


U<.utf!AL PAGE IS 

OF POOR QUALITY 


This subroutine generates the thrust THRUST and the fuel flow rate 
FOOT for some given altitude H, Mach number EMAKNO and EPR setting. First 
EPRMX, the maximum continuous EPR, is determined by table look-up for 
some given temperature Ta and altitude H, where 

T. - 1(1 + l§i(EHUa») 2 ) 2 - 273.15. 

Here* T is the temperature corresponding to altitude H, after temperature 
variation correction, and 

Y = 1.4, the ratio of specific heats. 

The EPR setting is limited to EPR <_ EPRMX for cruise and EPR < EPRMX - .1 
for climb or descent. 

Second, (FN/5 e ) is computed from a polynomial. Then, the thrust is 
computed as, 

THRST = 3(5 ) (FN/5 ). 

am e 

This is the thrust for the medium range tri-jet transport aircraft model. 
Here, 5^ is the pressure ratio 



Here, P is the atmospheric pressure corresponding to altitude H, and P q is 
the atmospheric pressure at sea level. A factor of three is used since 
there are three engines. 

Finally, the fuel flow rate FDOT is computed as: 

FDOT = 3*WFC * 5a *Kc 

where 

K - .00 223181 Ta + .9675897, 

c 

6 = 6 (1 + ^ ( EMAKNO) 2 ) y/y_1 . 

a am JL 


\ 

! 

1 


I 

.3 


1 



* 


] 

i 

1 

\ 
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Also, WFC Is the fuel-flow rate computed as a polynomial of EPR, where the 
coefficients of the polynomial depend on both altitude and Mach number. 


1 


Subroutines called: 


Co moons used: 


ATLOW 

DBLSRC 

P0LYE1 


CODE 

ERROR 

INPUT 

TRIJET 
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ENGEP3 


ORIGINAL PAGE fS 
OF. POOR QUALITY 


This subroutine generates the thrust THRUST and the fuel flow rate 
FDOT for some given altitude H, Mach number EM and EPR setting. First 
EPRMX, the maximum continuous EPR, is determined by table look-up , 

{Tables 6 and 7 in Block Data) for some given temperature Ta and altitude 
H, where 

Ta * T(1 (EM) 2 ) - 273.15. 

Here, T is the temperature corresponding to altitude H, after temperature 
variation correction, and 

Y * 1.4* the ratio of specific heats. 

The EPR setting is limited to EPR £ EPRMC. 

Second, (FN/5 ) is computed from a polynomial. Then, the thrust is 
e 

computed as, 

THRST =2(6 ) (FN /6 ) , 

am e 

where 6 is the pressure ratio 
am 



Here, P is the atmospheric pressure corresponding to altitude H, and P q is 
the atmospheric pressure at sea level. A factor of two is used since there 
are two engines. 

Finally, the fuel flow rate FDOT is computed. A polynomial is used 
to calculate WFC for a given EPR and altitude. At values of EPR < 1.6, 
there is also a correction for Mach number (Table 10 in Block Data). Then, 

FDOT = 2 *WFC * 5a *Kc , 

where 

K - -0022 T + 0.97 , 

c a 

6 a = 6 (1 + ^ (EM ) 2 ) y/y_1 . 

am 4 
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Subroutines called 


Counons used 



DfiLSRC 
POLY 2 
SGLSRC 


CODE 

ERROR 

INPUT 

TWINJT 



% 
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ENGIDL 


ENGIDL i'j called during descent to compute thrust and fuel flow rate 
for idle EPR. It does this through table look-up for the appropriate 
aircraft. 


Subroutines called: 


Commons used: 


DBLSRC ERROR 

SGLSRC TRUST 

TWINJT 


ESTCD 

This subroutine calls the appropriate routine to estimate cruise 
range before descent based on the particular aircraft model selected by 
the input variable IAC. Currently, two models are available, but logic 
is present to use up to four different aircraft. 

Subroutines called: Commons used: 

ESTCD1* None 

ESTCD2 

ESTCD3 

ESTCD4* 

* not included with program 


ESTCD2 

This subroutine estimates cruise range for the tri-jet, using a 
polynomial derived from experience with the trajectory computation. For 
the optimal profile, this polynomial gives a value of 115 nautical miles. 
For suboptimal profiles, it is less. 

Subroutines called: Commons used: 

P0LYE1 -BES-CRP 
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ESTCD3 


original. P A^T [y 
OE POOR QUALITY 


This subroutine estimates cruise range for the twin-jet. For an 
optimal flight, it estimates this range to be the difference between 
total range-to-go and climb plus the estimated descent range of 90 miles. 
For the suboptimal case, the cruise range is computed as a function of 
cruise energy and dl/d(energy) . 

Subroutines called: Commons used: 

3GLSRC CCDE 

CRUISE 

GRAPH 

DESCRP 

INPUT 


ESTDF 

This subroutine calls the appropriate routine to estimate the descent 
fuel based on the particular aircraft model selected by the input variable 
IAC. Currently, two models are available, but logic is present to use up 
to four different aircraft. 

Subroutines called: Commons used: 

ESTDF1* None 

ESTDF2 

ESTDF 3 

ESTDF4* 

* not included with program 


ESTDF2 — = = — — ■ 

This subroutine estimate descent fuel ior the tri-jet, using an 
emperically determined polynomial. 

Subroutines called: Commons used: 

None DESCRP 
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ESTDF3 


ESTDF3 estimates descent fuel for the twin-jet model as a function 
of the difference in energy between top of descent and landing. Different 
empirically derived functions are used depending on whether the single 
control or dual control case is being computed. 


Subroutines called: 


Commons used: 


POLY 2 CCDE 

CONST 

DESCRP 

INPUT 


ESTEP 

This subroutine computes the next energy step size during climb 
and descent. The nominal step is DENRGY, but this may be reduced when 
within 3000 feet of the estimated cruise energy. ESTEP also sets the 
quantities NEARCZ (within S00Q- feet of cruise) and IOPT (two-control 
flag only) . 

Subroutines called: Commons used: 


None 


CCDE 

DESCRP 

INPUT 


FBOUNP 


This function evaluates the drag force D in lb (using function DRAGC) 
if IDRAG is set equal to 1> and the absolute value of maximum thrust (TMAX) 
minus drag if IDRAG is equal to 2. To determine the maximum thrust, it 
calls the subroutine ENGEPR with thrust setting TMAX set equal to its 
maximum value* 

The output is: 

FBOUND « D in lb, if IDRAG = 1 

- | TMAX - D| , if IDRAG = 2 


Subroutines called: 


Commons used: 


DRAGC CCDE 

ENGEPR DESCRP 

ERROR 

INPUT 


FCLIMB 

This function evaluates the Hamiltonian for climb and descent for 
given values of airspeed VTAS, energy E, and EPR setting. The 
sB&titude is obtained by 

H = E - V 2 /2g. 

The atmospheric parameters and the aircraft ground speed V are found. 
The drag D is obtained by calling function DRAGC with Mach number and 
aircraft weight as inputs. The Hamiltonian is then evaluated by calling 
FTHRST . 


Subroutines called: 


Commons used: 


AT LOW 
DRAGC 
FTHRST 


CCDE 

CONST 

DESCRP 

ERROR 

INPUT 
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FCLMB6 


This function minimizes the Hamiltonian for climb. First SPLMT 
is called to calculate speed limits at the given altitude, and then 
the appropriate permissible region for minimization is computed. Then 
FCLMB6 calls FMIN and FCLIMB for minimization over speed. It also calls 
FMIN and FTHRST for minimization over EPR setting, if required for the 
two-control case (HT>I=1). 


Subroutines called: 


Conmons used: 


FCLIMB 

FDRAG 

FMIN 

FTHRST 

PI7-IMT 

SPLMT 


CCDE 

DESCRP 

ERROR 

INPUT 


FCOST 

This function evaluates the cost of flight per nautical mile, in 
dollars (cost of fuel and time) when IMFD is set equal to 0. It 
determines fuel flow rate in Ib/hr if IMFD is set equal to 1. To do 
so, FCOST calls first TRIM to obtain the trim condition for constant 
speed level flight at a given altitude H and the Mach number EM. 

Subroutines called: Commons used: 

TRIM CCDE 

WIND DESCRP 

ERROR 

INPUT 

WINDP 
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FDRAG 


This subroutine computes the drag D for some given true airspeed VTAS, 
energy E and aircraft weight W if IDRAG ■ 1. It computes T - D, if IDRAG t 1. 
Function DRAGC is used in computing D. The subroutine ENGEPR is called 
for computing T. 


Subroutines called: 


Commons used: 


ATLOW 

DRAGC 

ENGEPR 


CCDE 

DESCRP 

ERROR 

INPUT 


FDSCN6 

This function minimizes the Hamiltonian for descent. First SPLMT 
is called to calculate the speed limits at the given altitude, and then 
FDSCN6 calls FM1N and F CLIMB for minimization over speed. In the 
single-control case, this is all that is required. In the two-control 
case, FDSCN6 goes on to minimize over EPR setting as well. 


Subroutines called: 


Commons used: 


FCLIMB — 

FDRAG 

FMIN 

FTHRST 

P1LIMT 

SPLMT 


CCDE 

DESCRP 

ERROR 

INPUT 


FI AS 

FIAS returns Mach number as a function of indicated airspeed (in 
feet per second) and atmospheric pressure. 

Subroutines called: Commons used. 

None None 
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F1ASM 


FIASM returns indicated airspeed in knots as a function of Mach 
number and atmospheric pressure. 

Subroutines called: Commons used: 

None CONST 


FMIN 

This function minimizes the unimodal function F(X) so that x is 
within 1/144 of the range between input boundaries by a Fibonacci search, 
where 144 is the eleventh Fibonacci number. 


FMIN2 

This function minimizes the unimodal function F(X) so that x is 
within 1/34 of the input range by a Fibonacci search where 34 is the 
eighth Fibonacci number. 


koh* 


Thin (ami ton compute* the optimal nt vapeed toi minimum cruiae co» it 
it IMKP in not eoual to 0. It compute* opt ima l aitnpred tot minimum 1 %»«» I 
late It INKO in art equal I. Kivnt ! he minimum dtag 1>RUMN in computed 
bv calltug KMIN . Tlum. t ho maximum thiunt tn computed hv calling KNllKl'K, 
Thitd, the town and uppei hounds KA ami Kh. which dot too the i ange ot 
prtminnihlr Mach numhrm ate comput ed hv calling KMiN twice, timing 
KhOHNP tot evaluat ion, Klnallv the optimum Mach nurnhei Ol'TMAK and the 
minimum cost air computed hv (lie (unction KMIN, union Kt'OST tot evaluation. 


Suhvout tuen called; 

Arrow 

I'RAOO 

knoitk 

V'dOUNP 

tvosi 

KMIN 


Commons uned; 

tVPK 

PKSCRK 

KKKOK 

1NKKT 

WINOK 



V I NKS r 

Y\\\n t unci ton evaluates the Hamiltonian tot some given KKK netting, 
Mach numhet KM AON. tun* airspeed WAS, altitude UAl V. dtag 0 and weight W, 
Vhr values ot thiunt V and fuel t low tate KK ate obtained (tom calling 
the nuhioutiuc KNCITK. 



Suht out i tter. cal 1 ed ; 

Common* tmed 


KNutTK 

CCl'K 


KNuim 

tVNS V 


WlNP 

PKSCKK 

KKKOK 

g 


1NITT 

r- 


WINOK 


FULKST 


FULEST calls the appropriate routine to estimate the climb fuel 
use for the particular aircraft model selected by the input variable 
1AC* Currently, two models are available, but logic is present to use 
up to four different aircraft. 

Sub rout i n e s e ailed: Corarao ns used: 

FULEST 1* None 

FULKST 2 

KULST3 

KULKST4* 

* not included with program 


FI ULEST2 

This subroutine is used to estimate the climb fuel and the aircraft 
weight WCRIL! at the beginning of cruise for the tri-jet aircraft. WCKUZ 
is needed to interpolate the cruise table to obtain initial cruise 
values of energy, altitude, and cost. 

The routine uses cost of fuel FC, cost of time TC, initial weight WTO, 
and cruise, energy minus init ial energy to estimate fuel. 

Commons used: 

COW 
UESCRP 
KRROK 
INPUT 


Subroutines called: 
Wl.KFHV 


FULST3 


This subroutine is used to estimate the climb fuel and the air- 
craft weight WCRUZ at the beginning of cruise for the twin-jet aircraft. 
The routine uses cost of fuel FC, cost of time TC, initial weight WTO, 
and cruise energy minus initial energy to estimate fuel. 


Subroutines called: 

POLY 2 
WLEFHV 


Combo ns used: 

CODE 

DESCRP 

ERROR 

INPUT 


1 CLOCK 

This subroutine resolves TIME in seconds into hours IHR, minutes IMIN 
and seconds ISEC. For example, 4521 seconds is resolved into 1 hour 15 
minutes and 21 seconds. 


JTRUNC 


This subroutine finds the last point of a monotonically decreasing 
series . 
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LSQPOL 


This subroutine makes a least square polynomial fit for y as a function 
of x. The arrays x(0* y(’) contain corresponding points of the independent 
variable and the value of the function. The total number of points in the 
set is specified by N, the degree of polynomial to be fitted is specified 
by M - 1 (e.g. M ■ 3 for quadratic fit), and the coefficients of the poly- 
nomial are stored in the array B( • ) in the following order i 

y - B(3)x*x+ B(2) x + B (1) . 

LSQPOL uses the subroutine MATINV. Common used is ANE206 


MATINV 

This subroutine inverts a matrix A and stores the result in A. The 
dimension of the matrix is specified by M. The common used is ANE206. 


NICER 


This subroutine takes minimum and maximum values of the optimum 
trajectory variables and computes appropriate boundaries for the printer 
plots. 

NICER uses no other subroutines. 


PAGE 

This subroutine advances the printout to the top of the next page. 
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i i ^ 


Thin subroutine computes P. the percentage of variation trow 
The value of V Is computed by a local polynomial lit to the i> verau* K 
curves which have two points to begin (Rmin, 1.^4* or l.ty) and (Rmax, 1.014* 
or 1.04>). Subsequent Iteration tor synthesiflng the fixed range trajectory 
increases the number of points. The value of K and its corresponding 
percentage of variation are stored In the arrays RANOF and C respectively. 
The total numbei ot points at any Iteration is specified hv l PC. 

Subroutines called: Commons used: 

cm: 

DESCRV 
KRKOR 
CRAPH 
l NPUV 


FC1.FSV 
MAT \ N V 
WC1.ST 


o 


riClTR Ik, x, v, IS> 

This subroutine generates the printer plots and is called n -t * 
times itoi n points to be plotted} with k less than, equal to, or 
greater than rero, depending on the purpose of the call. 

First Call; k v 0. VICTVR stores one lx,v> point in the plot 
at ray and sets up the plot axes and scales. 


Next n Calls: k - 0, P1CTPR stores one (x,y} point In the plot 
array using the character designated by the para 
* meter IS. 


last Call: k > 0, PICTUK prints out the plot and writes the labels. 


Subrout l nes ca l l ed : 


Commons 


used : 


nn 


A 


N«'no 


r.RM'H 



PILIMT 


ORIGINAL PAGE IS 
OF POOR QUALITY 


This function evaluates the absolute value of T - D for some given 
EPR setting, altitude, Mach number and drag. It imposes a heavy penalty 
if the function is outside the boundary. Specifically, this penalty 
occurs for climb if it is outside the region where T < D. The penalty 
occurs for descent, if it is in the region where T > D. 

Subroutines called: Commons used: 

ENGEPR CCDE 

DESCRP 

ERROR 

INPUT 


P0LYE1 

This function evaluates the polynomial 

Y « b (1) + b (2) X + b(3)x 2 = . ..b(M)X™ _1 


P0LY2 


P0LY2 evaluates the polynomial 


Z - c u + c 12 x 2 + ... c lm x 2 


n-1 


c 21 x l + c 22 x l x 2 + •** c 2n x l x 2 


n-1 


+ .. 


t m-1 . m-1 , m-i. n-1 

+ c .x, + c 0 x- x* + ... + c x, x 0 

m2 1 m2 1 2 mn 1 2 
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PRETBL 


PRETBL prints out the cruise performance table for the lower cruise 
segment of a step climb flight. It uses the lower cruise table which 
was saved previously by subroutine STEP* \ 


Subroutines called: 


Commons used 

FI ASM 


CODE 

I CLOCK 


CRUISE 

SERCH1 


CONST 

WIND 


DESCRP 



ERROR 



GRAPH 


=r- 

INPUT 



TOA 



WINDP 


PRFTBL 

This subroutine writes the optimum climb and descent trajectory 
variables. It also writes the cruise performance table. In a step 
climb run* PRFTBL calls PRETBL to write the lower cruise segment and 
STEPUP to write the step climb trajectory. 


When the input quantity IGRAF is non- zero, this routine writes a 
dataset on Unit 11 for storage by the user. This dataset may be used 
for subsequent graphing. If IGRAF is greater than one, PRTPLT is called 
to produce printer plots. 


Subroutines called; 

Commons used 

CRUISR 

CODE 

CRUISX 

CRUISE 

FI ASM 

CONST 

1CL0CK 

DESCRP 

PAGE 

ERROR 

PRETBL 

GRAPH 

PRTPLT 

INPUT 

STEPUP 

WINDP 

WIND 
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PROFIL 


This subroutine controls the trajectory calculations. For a normal 
entry, it calls CR20P5 to calculate the cruise table and then calls PRSUM 
to print out the cruise summary table. PROFIL then goes through the 
following sequence. 

1. If this is a two-part trajectory, go to step 2. Otherwise, 
call CLIMB to generate the climb trajectory. 

2. If this is a constrained descent, go to step 9. 

3. Call DESCND to synthesize the descent profile, based on 
an estimated landing weight. 

4. Call VOPTRJ to estimate the cruise distance and fuel and to 
revise the landing weight. 

3. Call DESCND to generate the refined descent trajectory. 

6. Call VOPTRJ to generate the overall trajectory, including 
the total ground track distance covered by the trajectory. 

7. If this is a three-part, free altitude, non time-of-arrival run, 
call PCCMP3 to test this distance against the desired range R. 

If PCCMP5 returns with flag IRETRN, indicating that the trajec- 
tory has either been synthesized or cannot be, PROFIL returns to 
its calling program. If the trajectories to determine and 

^max ** ave not been synthesized, return to Step 1. If the re- 
quired distance is within the ir.terval (R , , R ), return to 

^ v min* max ’ 

Step 1. Otherwise if R is greater than RMAX, return to step 4, 

8. If PCCMP5 is. not called, PROFIL dimply returns. 

9. It a constrained descent is desired, call DESPC to calculate 
a new descent, and then return. 


Subroutines called: 


Commons used: 


CLIMB PCCMP5 
CRZ0P5 PRSUM 
DESCND VOPTRJ 
DESPC 


CCDE ERROR 

CONST INPUT 

CRUISE 
DESCRP 
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PR SUM 


PRSUM calculates and prints out the cruise summary table. It also 
calculates the suna«ed cruise time and distance tables. 


Subroutines called; 


Commons used: 


ATLOW 

PAGE 

WIND 


CCDE 

CRUISE 

DESCRP 

ERROR 

INPUT 

WINDP 


PRTFLT 

The subroutine calls subroutine PICTUR to generate printer plots 
of the optimum trajectory variables when the flag IGRAF is set greater 
than 1 in the subroutine PRFTBL. Labels to appear on plots are stored 
in array NOTES (8) which provides eight lines of ten characters each. 
The maximum and minimum values are calculated in subroutine PRFTBL and 
scaled in subroutine NICER. 

For IGRAF * 2, Mach no, flight path angle, altitude and 

fuel are plotted versus range. 

IGRAF ^ 2, all variables are plotted versus range and time. 

Subroutines called: Coimons used; 


NICER 

PICTUR 


GRAPH 

INPUT 


ORIGINAL PAGE IS 
OF QUALITY 


PRWT 



This subroutine prints the estimated conditions at top of climb. 


Subroutines called: 


Commons used: 


None 


CCDE 

CONST 

ERROR 


SERCHD 


The array TX(») is monoronically decreasing. This subroutine searches 
the index x such that 


Tx t+i 


9 


and returns both £• and pf where 


TX a ~ x 

TX^TX^i* 


SERCHI 

The array TX(0 la monotonically increasing. This subroutine searches 
the Index i such that 


TX t < X < TX f+1 , 


and returns both fc and pf where 
x - TX £ 

p ^ * tv : tv * 

TX i+ i 
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SGLSRC 


i 

iU 


: 5J. 










This function evaluates a single function F at the point A. This 
Is done by linear Interpolation to obtain A f s location in the array X 
and using the tabulated values of F(X) ♦ 

Subroutines called: Commons used: 

SERCHI None 


SPLMT 


SPLMT controls the speed limits for the aircraft during climb and 
descent. Different limits are applied during flight in the lowest 
10000 feet and the last 3000 feet. Limits also depend on the type of 
aircraft. 


Subroutines called: 


Commons used : 


ATLOW CCDE 

FT AS CONST 

SGLSRC DESCRP 

ERROR 

INPUT 





CP* 
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STEP 


STEP controls the logic for the step climb option. It is called 
(by either 0PTM56 or OPTTOA) after an initial optimum profile is computed, 
STEP then saves the initial cruise table and calls CRiI0P5 to calculate a 
new cruise table at an altitude 4000 feet higher. 

Function FMIN2 is then used to find an optimum trajectory consisting 
of a lower cruise segment, a step climb, and an optimum two-part cruise 
and descent segment, (See the description of STEPOPT for details of the 
function being minimized.) The new trajectory is only accepted if its 
total cost is less than the flight profile without the step dim;?. 


Subroutines called: 


Commons used: 


CRUISR 

CR20P5 

FMIN2 

PRSUM 

SERCHD 

STEPOPT 

STEPUP 


CCDE 

CRUISE 

DESCRP 

ERROR 

INPUT 

GRAPH 

TOA 


STEPEN 

Subroutine STEPEN computes the trajectory variables associated with 
a single energy step. It finds the altitude at the new energy level, 
updates the weather with ATLOW, finds the fuel flow from ENGEPR, and 
calculates time, distance, velocity, and flight path angle. 


Subroutines called: 


Commons used: 


ATLOW 

ENGEPR 

FIAS 

FIASM 

WIND 


CCDE 

CONST 

DESCRP 

ERROR 

INPUT 

WINDP 
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STEPOPT 


This function calculates a step climb trajectory. Input (from 
FMIN2) is the distance to be flown in the lower cruise table, STEPOPT 
then cruises this distance , calls STEPUP for a ramped Mach number 
climb at full power, and then calls PROFIL to compute an optimum two- 
part trajectory beginning at the top of climb* The final value of the 
function is the total cost of the four part trajectory consisting of 
lower cruise, step climb, upper cruise and descent. 


Subroutines called: Commons used: 


CRU1SR 

CCDE 

PROFIL 

CRUISE 

SERCHD 

DESCRP 

SERCHI 

ERROR 

STEPUP 

INPUT 

WLEFHV 

TOA 


STEPUP 


This subroutine computes a ramped Mach number climb at full power 
The climb is from altitudes HALT to HCRU in steps of DELTAH. __ 


Subroutines called: 


Commons used : 


ATLOW 

DRAGC 

ENGEPR 

F1ASM 

WIND 


CCD E 

CONST 

CRUISE 

DESCRP 

ERROR 

INPUT 

TOA 

WINDP 


TRACIT - _ 

In case of error, this subroutine provides a "walk back" through 
the subroutine calling hierarchy. If the subroutine is set up to 
recognize the computation or logic to be in error, then TRACIT can 
be used to find the source of the error* 
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T RIM 

This subroutine Ik used to compute the trim conditions (or 
medium range transport aircraft. This subroutine computes angle-of- 

at tack n and thrust T to keep the aircraft In trim for constant speed 
level flight , for a given altitude and for a given Mach number. 

With i the flight path angle, the equations of motion in the herisontal 
amt vertical directions are as follows: 

\ fdv/dt) • T cos u - P - W sin \ 

W 

v(d)/dt) * T sin a + L - VI cos •) 

For a t r imme J cond i t Ion : 

(dv/dt) - (awao - 0. 

Tho two equations aro combined to eliminate thrust to give the equation: 
(Wcosy - 1 ) cosvt - (Ositut + yja.iny)alnrt * 0, 

This equation is solved hv iterating with angle-ot -at tack , a. 

Oiuo the alrcratt Is trimmed, the thrust Is solved from the drag 
T * l P + Ws ill) ) / coso . 

This required thrust Is matched by iterating on values ot power setting 
(EFR) and calling subt online ENUKTR. Once the correct power setting is 
determined, the engine fuel How is also known* 

Subroutines called: Commons used: 

CUKAi: COOL 

cum DEscnr 

ENU1TR ERROR 

INPUT 


vornu 


This sub t out t no {\ ) eomputos t ho ollmh ft mo 1 KCtMB, ollmh t Imo TCClJMtt 
;* i imti dlsiauoo IH'IWH, dosoonl tuol KOOWN. dosoont t imo THOVIN, mid drwcont 
distauoo UI10WN; (?) ovulso tuol KCKUl.H, ovulso distauoo IH’KUZ, mid oiulso 
time TORU?.; ( D tho final o rut sc' wo(£hi WCKU? ; (*0 I ho oiulso tuol 
ottiotouov KKCRIU*. (Ih/u.iui. ) and ovoiall tuol otiloionov RKKl'N?. (Ib/n.mM) 
and (<») t ho hnultMji wot Kit t Wl.NDC. It also Ronotat os tho tuol usod hi 
lh. the' distauoo oovovod in mill's, tho time' durst ion tu hi : mill: soo. 
tho toial oust in $, ami tho oust poi nunlioal mtlo tot ollrab* ovulso, 
dosoont • and tho ovoiall tva|ootnrv. 

It VOriK.l is oallod dui tiiK stop ollmh opt tmt /at Ion , it iucludos 
tho lowoi o tu t so aoRinout and stop ollmh timo, tuol and dlstanoo in 
t hoso totals. VOPTR.1 also prints out tho initial orutso woiRht , final 
ot ui so wotRht. amt tholi ooi i ospondiiiR tvuo altspood, ovulso oust, 
oqutvalont aitspood, oiutso ouoirv, Rt ound spood , altitudo and Maoh numlvr 
Ktnallv, it pt tuts out tho ovulso tuol olfioionov, tho ovoiall tuol 
ottiotouov. and tho landing, woiRht » 


Suhvoutinos oallod; Commons usod: 


ATI. OK 

STKCtll 

caiK 

KKKOR 

aunsx 

W l NO 

CONST 

OKAHl 

l‘l ASM 

W1 KK»V 

CRIUSI* 

INPUT 

»’Ai:v 

WKMTl 

I'KSORP 

viNor 


Ml 


WATKST 

This subroutine is used to estimate the cruise range, iuel burned 
during descent (FIXWN) and the aircraft landing weight Wl.Nlx:. For tree 
cruise altitude optimization using V and n, (1VF1 - l) the cruise range 
is zero. For optimisation using V onlv, the cruise distance is estimated 
using KSTOU, For fixed cruise altitude, the cruise distance is assumed 
Co be the range minus the climb distance minus ^0. CKtUSX is called 
to bring i he flight parameters to this range lest tmated top of descent 1 , 
Then WI.FFHV is called to estimate cruise fuel rate and ground speed. 

For all cases, KSTOF is then called to compute the descent fuel. 

The landing weight is then calculated from climb, cruise, and descent 
fuel vi . jo . 


Subroutines called: 

Commons i 

CRIHSX 

tXOF 

FSTCP 

CONST 

F.STl'F 

PFSCKT 

Wl.FFHV 

FRROK 


CKAHI 


XWVV 


TO A 


W [01. ST 

This subroutine is vised to compute climb, ciuisc, and descent 
segments toi the fixed thrust case ilVPl - 0^ where the value of 
cruise cost is very close to the optimum value. This toutino is called 
from IHWb. 


Subtout lues called: Commons vised: 


Cl l Mb 

lYPK 

PFSCNP 

CONS )' 

pf.sfc 

nvst'Ki 

vornu UU entvv point V01M 

KKKOK 


GRAPH 


INPUT 


u: 


WIND 


This subroutine computes the wind velocity os a function of altitude. 
This is combined with the aircraft velocity with respect to the air mass 
to compute ground velocity. Inputs to this program are H* the altitude 
in feet; PSIG* the aircraft heading in degrees; VTAS the aircraft air 
speed; GAMMR* the angle of attack; Wl* and PSIW arrays. The outputs 
from this program are WA, the wind speed* and VG, the aircraft ground 
speed. 

Subroutines called: Commons used: 

SERCHI DESCRP 

ERROR 

INPUT 

WIND? 


WINDIN 


This subroutine reads in the wind profile (the magnitude and direction 
of wind as a function of altitude). The wind magnitude input is in knots 
and the program converts it to ft/sec and stores it in the VW array. The 
wind direction is stored in PSIW in degrees. The input represents the 
direction the wind is coming from. The program adds 180° to this value 
to obtain the vector direction. The altitudes corresponding to these 
wind magnitudes and directions are stored in array HWIND. The wind may 
be input as a single profile valid over the entire flight* or as separate 
climb* cruise and descent profiles. In the case of a step climb the cruise 
profile is used for lower cruise* step climb, and upper cruise segments. 

Subroutines called: Commons used: 

None INPUT 

V1INDP 


U 3 


WLEFHV 


The purpose of this subroutine is to interpolate cruise table results 
to obtain cruise cost ELAMBS, cruise energy ECRUZ , cruise fuel flow rate 
FCRUZ, cruise altitude HCRUZ, and c.ruise airspeed and ground speed VCKTAS 
and VGKNT as functions of the input cruise weight WCRUZ. 


Subroutines called: 

Commons used 

CTABLE 

CCDE 

JTRUNC 

CONST 

SERCHD 

CRUISE 

WIND 

DESCRP 


ERROR 


INPUT 


WINDP 


WRITE1 


This subroutine (a) computes the climb, cruise, descent, or overall 
cost for flying the specified segment of the trajectory, where the cost 
Is the sum of fuel cost and time cost, (b) computes the cost per nautical 
mile to fly, (c) resolves the time given in seconds into hours, minute 
and seconds, and (d) writes out the description of the segment (i,e. climb, 
cruise, descent, total), the fuel used, the distance traversed, the time 
duration (in hours, minutes, seconds), the cost for the segment, and the 
cost per nautical mile to fly. 


Commons used: 


Subroutines called: 
ICLOCK 


CONST 

INPUT 
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